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ABSTRACT 

Collisionless entrainment of the surrounding matter imports the relativistic baryon component in the Gamma- 
Ray Burst (GRB) fireball frame. We show that half the fireball energy can be transferred from radiation to the 
comoving hot motions of baryons under the photosphere. The yet baryon-poor fireball can reexpand to a 
very high Lorentz factor (VHLF) T ~ 10^-10^ by its own relativistic collisionless pressure beyond the photo- 
sphere (so-called collisionless bulk acceleration), leading to internal and external shocks. A simple synchrotron 
emission from the VHLF internal shocks produces (i) the extra power-law spectral component with variability 
observed in the Fermi Ge V bursts, up to the Te V range for the future Cherenkov Telescope Array (CTA), (ii) the 
GeV onset delay with a weak luminosity dependence fdeiay ^ L~^^^, and (iii) the spectral break of GRB 090926 
by the synchrotron cooling break or the maximum synchrotron cutoff limited by the dynamical time, not by 
the creation cutoff. The relativistic baryon component could also heat the photospheric thermal photons 
into the main GRB Band spectrum via pp, p-f (Bethe-Heitler and photomeson), and Coulomb thermalization 
processes. In this hot photosphere-internal-external shock model, we can predict the anticorrelation of ^TeV 
neutrinos and GeV 7-rays, which may be detectable using IceCube. The spectral peak and luminosity (Yone- 
toku) relation is also reproduced if the progenitor stars are nearly identical. We also discuss the steep/shallow 
decay of early X-ray afterglows and short CjRBs. 

Subject headings: gamma rays; bursts — gamma rays; theory — radiation mechanism; non-thermal 



1. INTRODUCTION 

The cosmological Gamma-Ray Bursts (GRBs) are the 
most luminous objects in the Universe. Although the 
physical processes at the central engine are far from un- 
derstood, the fireball model is generally accepted as the 
paradigm for produ cing the relativistic outflows and high- 
energy em ission Cavallo & Rees' '1978'; 'PaczvhskJ '1986t 
Goodman Il986t ISh emi & Piran 1990; Meszaros & Reesl 
1993aiibt iGrimsrud & Wassermani il998t iMeszaros & Reesl 



2000: 



However, the actual emission mechanism of prompt GRBs 
is still debated, lacking a consistent picture. The main prob- 
lem is the high efficiency (> 50%) of the GRB prompt emis- 
sion, defined by the GRB ener gy divide d by the total energy 
including the afterglow energy jZhang e t al. 2006; lok a et all 
I2OO6). With the use of the internal shock model (the leading 
model for the prompt emission), it is difficult to achieve a high 
efficiency wi thout a large disp ersion in the Lorentz factor of 
the outflow s ('Kobav ashietaI]ll997tlKobayashi & Sarill2001l; 
iBeloborodo v 2000). Even if it is achieved, a large dispersion 
in the Lorentz factor makes it difficult to realize the observed 
spe ctral relation s (Zhang & Meszaros 2002), such as the Am- 
ati dAmati. Frontera & Guidorzi 2009; Amati 2010) a nd Yo- 
netoku relations dY^netokuet aL.2004: .Kodama et al J 120081; 
iNava et aUlMoh . 
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where Cpeak is the peak energy of the observed broken power- 
law spectrum (so-called the Band spectrum) and L is the ap- 
parent isotropic luminosity of the prompt emission observed 
within an angle < l/F of a jet axis. S uch a correlation 
is also satisfied with in individual pulses (lOhno et al.l 120081: 
iGhirlanda et al.ll2()09b . Another problem of the internal shock 
synchrotron model is that the low-energy spectral slope be- 



comes steeper than that observed, owing to the intrinsic syn- 
chrotron spectru m and, e ven worse, owing to the fast electron 
cooling (Ghisel hni et al.l l2000a; Me szaros & R ees 2000a). 

These difficulties of the internal shock m odels lead to the 
reexamination of the original fireball model dPaczv nski 198^ 
iGoodman 1986). in which photons are released as photo- 

fheric emission when the fireball becomes optically thin 
hompson 1994; Rees & Meszaros 2004; Thompson et ap 
2007; " Meszaros & Rees. 20 00a; De rishev et al. 1999[ 
Meszaros et al. 200 21; lPe'er& Waxman r|2'004t Thompson! 
2006; Pe'eretal. 2006^ 'Rvde 2005; 'ciannios & Spruil 
2007; Pe'eretal. 2007; Giannios 2008; lokaetal. 2007r 
Ghisellini et al. 2007; Rvde & Peleij I2009t (Beloborodo^ 
,2009; Lazzati & Begelman 2010)- The original problem 
is alleviat ed by introducing the d is sipation under the pho - 
tosphere (i Rees & Meszaros! I2004t iTTiompson et al.l 120071) . 
which can bring the thermal peak into the observed range in 
Eq. ([T]i. The photosphere model can naturally achieve the 
high efficiency and the hard low-energy spectrum. The only 
flaw is that the spectrum tends to be thermal without nonther- 
mal tails observed in GRBs, although a substantial fraction 
(~ 30%) of long GRBs may have thermal peaks (RvcS 
l2005l;lRvde & Pe'ei!!2009!) . The nonthermal tails could arise 
from Comptonization of the thermal photons by electrons 
and positrons 
reconnection 



(e^), heated at diss ii 



ation, such as m agnetic 
IGianniosI "2008), 



pation , 
pOTl 



Gi an nios & Sprui 
neutron collisions (Derish ev et al.lll 999; Beloborodov 2009), 
or re peated shocks (loka et al.l 12007; .Lazzati & Begelman! 
120101) . However, it seems questionable that a subdominant 
or a different component rather than the thermal component 
can supply the nonthermal energy that is accidentally com- 
parable to the dominant thermal energy. Nevertheless, the 
photosphere model has an advantage that the peak energy 
Epeak IS Stabilized as it is fixed by the temperature of the 
photosphere, regardless of the dissipation mechanism. The 
spectral relations in Eq. ^ are more easily reproduced 
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dRees & MeszarosI l2004t iThompson etaLl l2007h than the 
other attempts t o solve the emission mechanism, such as th e 
jitter radiation dMedvede v & Loeb" 19991 'Medvedev"200q^ 
Klein- Nishina effect jP erishev et al. 2000; Bosniak et 

f)9; Wan2 et al. 2010), synchr otron self-Compton (SSC) 
inaitescu & Meszaros 2003l; IStern & Poutanen 2004, 
bulk Compto n (Lazzati et al. 2000; Ghisellini et al. 2000fi 
Lazzati et al.l |2004|). decaying magnetic field CRossi & ReesI 
20031; iPe'er & Zhand 120061). a nd q u asi-thermal Comp- 
tonization (iGhisellini & Celottil 119991; lAsano & Terasawal 
12009). 

Recently, the Fermi satellite, launched on 11 June 2008 
with the GBM (8 keV - 40 MeV) and LAT (- 20 MeV - 
300 GeV) detectors, has been used to observe ^GeV 7-rays 
from GRBs, providing interesting clues to the emission mech- 
anism (Abdo et al. 2009a b.c d, 2010a b.c). The GeV events 
are increasing more than sixfold from the era of EGRET that 
detected an 18 GeV photon 90 min after the burst in GRB 
940217 (lHurlevetal.(ll994l) and a rising late Ge V spectral 
component in GRB 941017 jGonzalez et alj|2003l) . The main 
features of the Fermi bursts are summarized as follows: 

(1) In some Fermi/LAT bursts, the Lorentz factor of the 
outflows is constrained in the relatively high range 
r > 10^ so that the high-energy p hotons can avoid th e 
annih ilation by pair creation jAbdo et al.ll2009b llal 
l20T0d) . 

(2) Fermi found an additional spectral component at > 10 
MeV with comparable energy to that of the main Band 
component, at least, in short GRB 090510 (Abdo et al.' 
2010c.) and long GRB 090902B ( Abdo et al. 2009a; 
Ryde et ani2010l) . This extra component is fitted by a 
single power-law that slightly rises in i/F^, oc J/" '-!/"'* 
and often extends below < 20 keV over > 7 energy dig- 
its. (However, note that no other experiments have con- 
firmed the low-energy extension of the power-law com- 
ponent.) 

(3) The high-energy (> 100 MeV) emission lasts longer 
than the MeV emission in most LAT GRBs (lAbdo et alJ 
|2009d b, 2010a, 2009a, 2010c). The well-observed ex- 
tended emission shows a temporal power-law decay up 
to > 10^ s. 

(4) The high-energy emission sometimes sh ows a large 
amplit ude variability on short timescales dAbdo et al.l 
l2009ah . 

(5) The high-energy emission is delayed behind the on- 
set of the MeV emission i n almo st all LAT GRBs 
(lAbdo et a l. 2009d,bt 1201 Oai 120091 120103) . The de- 
lay time in the rest frame is fjeiay ~ 1 s for long GRBs 
and - 0.1 s for short bursts, GRB 081024B and GRB 
090510. 

(6) The prompt emission spectrum of GRB 090902B has 
a quasi-blackbody component, which is consistent 
with the photospheric emission (lAbdo et al.1 l2009at 
iRyde eTaL-^OlO) . 

The Fermi discoveries excite the theoretical reconsid- 
erations of the GeV emission, which may be classified 
as (i) the external shock models with synchrotron emis- 
sion from adiabatic shocks (.Kumar & Barniol Duraa2009a.b; 



Barniol Duran & Kumadl2010HCorsi et al.l2 009allbl;lGao et alj 
2009t [De Pasquale et al. 2010; Pandey et al ll2010l) and ra ^ 
diativ e shocks '(.G hirlanda et al. 2010; Ghisellini et al."2010l), 
SSC (IZou et al] 2009; Wang et al. 2010; Corsi et al. 2009^ 
lNeamusll201 (fr and external Compton (Murase et al. 2010), 
and (ii) the internal shock models with synchrot ron emis- 
sion (Wang et al. 2009; Fa nll2009b. SSC (Li 2010a: IZou et alj 
l2009i:,Abdo et al. 200 9a; Corsi et al."2009a b), hadronic emis- 
sion (lAsano et al.f|20 09a; Razzaque et al. 2009; Asano et al] 
l20091j) . and external Compton of cocoon (Toma et alji2009l) 
or photospheric emission (Toma et al. 2010). The extended 
GeV emission most likely has the external shock origin. The 
emission mechanism could be other than synchrotron since 
the maximum synchrotron cutoff terminates the late (> 100 
s) emission o f > 10 GeV photons from external shock syn- 
chrotron (O I2010bt Piran&Nakar 2010). On the other 
hand, the external shocks cannot produce t he observed large 
amplitude variabilit y on short timescales (lloka et alJ l2005t 
ISari & Piranlll99'7bl) . so that an additional origin, probably 
the internal shock emission, is also required. However, with 
the use of internal shock models, it is difficult to explain 
the hard extra component that extends to the low-energy ex- 
cess. Obviously, synchrotron emission only cannot produce 
the extra component in addition to the Band component. 
The SSC emission usually peaks at higher energy than syn- 
chrotron without extending to the low-energy excess. Al- 
though hadronic models can make a low-energy excess via 
direct and cascade radiation (e.g., synchrotron emission by 
secondary pairs at low energies), the proton injection isotropic 
luminosity should be larger than 10^^ e rg/s, posing a challenge 
for these model s (| Asano et ani2()09bl) . The external Comp- 
ton of cocoon dlbma et alj 1200 9") or photospheric emission 
(lTomaetal.ll2010ir seems viable, but would need a fine tun- 
ing to smoothly connect the high- and low-energy excesses 
that have different emission origins in these models. There- 
fore, the prompt emission, not only the main Band component 
but also the extra component, remains a mystery. 

In this paper, we revisit the dissipative photosphere model 
in light of the Fermi results, scrutinizing the dissipative pro- 
cesses that reproduce the spectral peak and luminosity (Yo- 
netoku) relation in Eq. dTJ. We suggest that the dissipation is 
caus ed by the surrounding matter that decelerates the fireball 
(see ^3.1l l. We show that the matter is loaded as the relativistic 
baryon component into the fireball, with a significant amount 
of energy received, comparable to that of the radiation com- 
ponent, using the energy and momentum conservation (see 
^3.2b . This hitherto missing component can alter the fireball 
dynamics and spectra in a favorable way to solve the prompt 
emission. Firstly, if not fully thermalized as expected in the 
baryon-poor fireball, the relativistic baryon component can re- 
expand to a very high Lorentz factor (VHLF) F - 10^-10^ by 
its own relativistic collisionless pressure, via a so-called "col- 
lisionless bulk acceleration" mechanism (see ^for a short 
summary). The subsequent VHLF internal shocks can explain 
the extra high-energy component with variability in the Fermi 
bursts by a single emission mechanism of synchrotron (see 
®. Secondly, the relativistic baryon component could also 
operate as a heating source for to Comptonize the pho- 
tospheric thermal photons into the observed Band spectrum 
via pp, p^ (Bethe-Heitler and photomeson), and Coulomb 
thermalization processes. Without fine tuning, the relativistic 
baryon component has the right amount of energy (compara- 
ble to the thermal energy) to make the nonthermal tails (see 
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§ 33.21 and l5. 21 ). In a sense, we consider a "hot photosphere" 
that is only partially thermalized with relativistic relic parti- 
cles. Our picture fall s into th e photosphere-internal-external 
shock scenario (To ma et"an i2010). in which the main Band 
emission comes from the photosphere and the extra compo- 
nents come from the internal and external shocks. 

This paper is organized as follows. In ^ we first sum- 
marize the basic idea for making the VHLF fireballs. In ^ 
we go into the fireball dynamics, expanding t he id ea of the 
collisionless bulk acceleration to a VHLF. In 33.11 we recall 
that the observed Yonetoku relation in Eq. ([U strongly sug- 
gests the fireball dissipation under the p hotosphere, probably 
caused by the baryon loading. In 33.21 we use a simple two- 
body collision to describe the dissipation, properly taking into 
account the relativistic hot motions of baryons before ther- 
malization. We show that the relativistic baryon component 
naturally achieves comparable energy to the radiation. T hen, 
after giving the photospheric and pionosphe ric ra dii in 33.31 
we derive the final coasting Lorentz factor in 33.41 In 33.51 we 
examine the thermalization processes via p7 (Bethe-Heitler 
and photomeson). Coulomb, and plasma interactions in addi- 
tion to pp collisions. We devote 33.61 to consistency checks 
with previous works, and 33.7l to remarks on the connections 
between the collisionless bulk acceleration and the magnetic 
acceleration in making the VHLF fireballs. 

Secondly, in ^ we apply the VHLF fireballs to the extra 
hig h-en ergy component with variability in the Fermi bursts. 
In 34.11 we firstly argue the GRB emission site in fav or of 
the photosphere-internal-extemal shock scenario. In 34.21 
we calculate the internal shock synchrotron spectrum in the 
VHLF models, which is found to be consistent with the ob- 
servations, and also a nice target for the future Cherenkov 
Telescope Array (CTA), since the creation cutoff goes be- 
yond the TeV range in the VHLF models. In 34.31 we discuss 
the possible origins of the spectral break around ^ 1 .4 Ge V 
observed in the extra component of GRB 090926, which was 
suggested as the creation cutoff for F ^ 600. In the VHLF 
models, the spectral break could be the synchrotron cooling 
break for T - 10^ or the maximum synchrotron cutoff lim- 
ited by the dynamical time, for F 10^. In 34.41 we sug- 
gest that the GeV onset delay and its weak dependence on the 
burst parameters can be naturally derived if the baryon load- 
ing at the dissipation is rich shortly after the jet break out of 
the progenitor star. 

Finally, i n ^ we discuss the future predictions and open 
issues. In 35.11 we predict the anticorrelation of ^TeV neu- 
trinos and GeV 7-rays, which might be detectable using Ice- 
Cube. In 35.21 we suggest that the relativistic baryon compo- 
nent could transform the photospheric thermal photons into 
the observed Band spectrum. In 35.31 we go back to the ori- 
gin of the Yonetoku relation in Eq. ([TJ, which implies that the 
baryon loading rate, i.e, th e env ironmental condition, is nearly 
identical to any bursts. In 35.41 we discuss possible configura- 
tions of the baryon loading. We also discuss the model impli- 
cations for the steep/shallow decay of early X-ray afterglows 
in 35.61 and short GRBs in 35.71 We use the cgs units with 
kB = h=\ and = Q/W, and the standard cosmological pa- 
rameters, Q,„i = 0.3 and VLjy = 0.7, unless otherwise stated. 

2. BASIC IDEA FOR VERY HIGH LORENTZ FACTOR (VHLF) 

In this section, we briefly summarize the essence of the 
VHLF fireball formation, and clarify the connections between 
following sections. 

Let us quickly recall the conventional fireball model 



(iPaczviiskil [ 1986^ lGoodman| [19861 iShemi & Piraril fT990t 
iMeszaros & Rees 2000a). We initially consider a fireball that 
is compact and radiation-dominated. Since it is optically 
thick, the fireball expands by its own radiation pressure. The 
Lorentz factor grows as F cx r, and is saturated at a constant F^ 
after almost all the radiation energy is converted into kinetic 
energy. The coasting Lorentz factor Tc is higher for lower 
baryon loads, whereas it has an upper limit for sufficiently low 
baryon loads because the fireball becomes optically thin in the 
accelerating phase before converting all the radiation energy 
into kinetic energy (i.e., before the saturation). The maximum 
value of the coasting Lorentz factor is about Tc < ?/* ^ lO"* for 
typical parameters [see Eq. ( |23] | and 33.31 . 

Our idea is that even the baryon-poor fireball can accelerate 
to the saturating (very high) Lorentz factor, if the radiation 
energy is transferred to the other relativistic component that 
is not radiated away from the fireball. In our case, this is the 
relativistic baryon component (and the magnetic component 
resulting from it). We initially consider an almost baryon-free 
fireball (e.g., a leptonic fireball). A small amount of baryon 
is loaded into the fireball as the relativistic baryon compo- 
nent from the surrounding matter near the progenitor star or 
the preceding baryon-rich ejecta. The loaded baryon decel- 
erates the fireball, receiving energy from the radiation in a 
collisionless way. We show that half the radiation energy can 
be transferred to the relativistic baryon component at the fire- 
ball dissipation using the energy and momentum conservation 
in Eqs. ( flTT i and ( fT2b . 

The fireball dissipation can occur under the photosphere 
(before the fireball becomes optically thin) so that the dis- 
sipated energy is trapped by the fireball. The dissipation is 
even demanded by the GRB observations if we identify the 
spectral peak energy £peak with the fireball te mperature (see 
33. H and Refs. lRees & Mes zaros (2004) and Tho mpson et alj 
(2007)). The observed temperature Fr' drops at the dissi- 
pation because the fireball radius at the base of the flow r/, 
effectively increases, and this enables us t o ac cou nt fo r the 
observed spectral relations in Eq. ([T]) (see § 33.1l and l5.3l ). The 
thermalization of radiation and is much faster than that of 
protons, keeping the comoving leptonic temperature at a non- 
relativistic value. 

In previous studies, we usually (implicitly) assume that the 
baryon (proton) component of the fireball is also completely 
thermalized at the dissipation. This is the case (A) in Figs.[T| 
and 121 where we schematically show the fireball evolution 
with dissipation. The radiation-dominated fireball expands 
outward with F oc r, and is dissipated, e.g., via shocks with 
baryon (protons) initially at rest. Here, if we look closely into 
the dissipation, the shock is usually collisionless at least in the 
early stage. The protons are trapped by the fireball via mag- 
netic fields, which could be generated by the plasma instabili- 
ties or could be advected from the central engine. In any case, 
the velocity of protons is not changed so much in the frame 
of the shocked region. Then, the isotropized protons have a 
random Lorentz factor of 

7; ~r,„ (^10^-10^), (2) 

in the shocked region because the protons drive into the 
shocked region with the Lorentz factor that is about the bulk 
Lorentz factor of the fireball after the merger F,,,. In this pro- 
cess, the fireball energy is transferred to protons, leading to 
an almost equipartition of the (comoving) energy density be- 
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Fig. 1 . — Schematic of a fireball evolution with dissipation. The radiation- 
dominated fireball is decelerated to a Lorentz factor F,,, ~ ^yTj under the pho- 
tosphere, via shock with baryon (protons) initially at rest. The comoving 
energy density of protons can be boosted to a value comparable to the radia- 
tion. Up ~U!y, where almost all the proton energy is carried by the relativistic 
hot motions of protons with a random Lorentz factor of 7^ ~ F„, . (A) If ther- 
malization is effective, almost all the proton energy immediately dissipates 
into radiation, leading back to a standard radiation-dominated fireball. The 
maximum Lorentz factor is less than r;» ~ 10^ as in the conventional case. 
(B) If thermalization is not effective, the relativistic collisionless motions of 
protons continue to push the fireball with F oc r up to a (saturating) VHLF, 
~ Yp^m ~ r; ~ 10^-10*, even beyond the photosphere. The kinetic en- 
ergy remains comparable to the radiation energy since the energy density of 
relativistic protons behaves like radiation, Up ~ C/^ oc r"*. 

tween protons and radiation [see ^3. 21 and Eq. (fTSlll. 

u'p ^ f/;, (3) 

according to the energy and momentum conservation in 
Eqs. ( fTTT i and (fTZb . Thus, the kinetic luminosity temporar- 
ily equals the radiation luminosity, L^iA) ^ L^{A), in Fig.|2] 
The energy equipartition also renders the bulk Lorentz factor 
after the merger as 



r„ 



(4) 



where 77 is the dimensionless entropy (the radiation to baryon 
ratio) of the fireball after the merger in Eq. fT% . How- 
ever, almost all the proton energy is carried by the relativis- 
tic hot motions of protons, not by the rest mass energy, i.e., 
^prest ^ ^p- The subsequent thermalization of protons via 
pp, pj (Bethe- Heit ler and photomeson), and Coulomb in- 
teractions (see 33.51 ) reduces the proton energy and pressure 
considerably into radiation [For example, the protons are ef- 
fectively thermalized via pp collisions under the pionosphere 
where the optical depth to pp collisions is larger 



r < r, 



pp- 



than unity Tpp > 1 (see 33.3b 1. Then, the fireball evolution 
is essentially similar to that without shocks since the dissi- 
pated energy is trapped in the system under the photosphere. 
Therefore, the coasting Lorentz factor has a conventional up- 
per limit of 



10^ 



(5) 



as the radiation escapes from the photosphere before transfer- 
ring its energy to the baryon kinetic energy [see Eq. ( l23T l and 
El. 

The evolution is totally different if the thermalization is not 
completed. This is the case (B) in Figs. [T] and |2] The evo- 



> 5r 

yj 9- 

O 



10° 
10^ 



10' r 



N lO'' 



U/B) 



N il i I 



Lk(A) i: 



4- 



4- 



Lk(B) 



4- 



4- 



Intemal Energy 
is thermaiized 
into Radiation 



(B) fshock* Tppi 




10" t- 
10^ 



10" 



Fig. 2. — Schematic evolution of physical quantities as a function of the 
fireball radius r. The radiation-dominated fireball is decelerated via shock 
with baryon (protons) initially at rest, boosting the proton energy to a value 
comparable to the radiation. (A) (dashed lines) If the thermalization is effec- 
tive, e.g., under the pionosphere r < rpp with the pp coUisional optical depth 
Tpp > 1 , almost all the proton kinetic luminosity Li^(A) immediately dissipates 
into the radiation luminosity Lj(A). The Lorentz factor grows under the ra- 
diation pressure, so that it coasts at the photosphere r = Tpj, with F^ < lO'. 
(B) (solid lines) If the protons are not thermalized, the kinetic luminosity re- 
mains comparable to the radiation luminosity, Li-(B) ~ L-y(B), since almost 
all the proton energy is canied by the relativistic random motions of protons, 
not by the rest mass energy, Z-t.restlB) Lk{B). The Lorentz factor grows 
under the proton collisionless pressure up to a VHLF, F^ ~ F,^, ~ lO-^-lO', 
even beyond the photosphere r > rph. In both cases (A) and (B), the observed 
temperature TT' drops by ~ (rt/ro)^^' times because the fireball radius at 
the base of the flow r/, effectively increases. 

lution is almost the same as the previous case (A) before the 
proton thermalization. However, if the relativistic hot motions 
of protons are not thermalized (e.g., pp collisions are not ef- 
fective with Tpp < 1), the relativistic collisionless motions of 
protons reexpand the fireball with F cx r, acting like radiation 
pressure. Since the random motions are converted into the 
bulk motion, the final bulk Lorentz factor is the bulk Lorentz 
factor after the merger r„, multiplied by the random Lorentz 
factor 7' ^ r,„ as 



77 



10"*- 10" 



(6) 



which is a (saturating) VHLF, much larger than the conven- 
tional upper limit 77* ~ 10^ in Eqs. (|5]l and ( l23l l. It is remark- 
able that such a "collisionless bulk acceleration" continues be- 
yond the photosphere, r > rph, i.e., even after radiation loses 
contact with matter, in contrast with the thermalization case 
(A). The kinetic luminosity also remains comparable to the 
radiation luminosity. 



(7) 



since the energy density of relativistic protons behaves like 
radiation. Up ^ U!y (x r~^. 

The energy source of the collisionless bulk acceleration is 
the initial radiation energy. The radiation energy is transferred 
to the relativistic baryon component in the nonradiative form 
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of the isotropic random hot motions, which is later converted 
into the bulk kinetic energy. Even if the relativistic baryon 
energy dissipates into the magnetic field, a similar accelera- 
tion continues via the magnetic pressure (see § ^3.5l and l3.7l i. 
However, the acceleration mechanism does not work if the 
baryon loading occurs after the coasting since adding mass 
just reduces the bulk Lorentz factorQ 

According to the above considerations, it is physically rea- 
sonable to define a VHLF as a Lorentz factor larger than the 
conventional maximum value 77, ^ 10^ in Eqs. (|5]l and (|23] |. A 
VHLF could open a new paradigm for interpreting the GRB 
properties, removing a theoretical bias to the lower Lorentz 
factors. In ^ we apply the VHLF fireballs to the inter- 
nal shock synchrotron model for reproducing the keV-GeV 
power-l aw s pectrum with high time variability detected by 
Fermi (EH, the spectral break in GRB 090926 by the syn- 
chrotron cooling break or the maximum synchrotron cutoff 
limit ed by the dynamical time 04.3t . and the GeV onset de- 
lay ( ^4.41 ). In H we predict an anticorrelation between GeV 
7-rays and TeV neutrinos, and also suggest that the relativistic 
baryon component could heat the photospheric thermal pho- 
tons into the observed Band spectrum. 

Note that each element of our idea to create VHLF fire- 
balls is not completely new. The acceleration, which converts 
the internal energy back into the kinetic energy, was previ- 
ously discussed in the context of the int ernal shock efficiency 
dKobavashi & SarilllOOlt lKumari[T999h . although their fire- 
balls cannot reach a VHLF since shocks occur in the coast- 
ing phase, not in the accelerating phase. The fireball dissi- 
pation under the photosphere was also di scussed in the photo- 
sphere mode l (Rees & Meszaros 2004; T hompson et al.l2007l : 
iGhiselUnfet al.i i2007i: Joka et al.i i2007D . The pionosphere 
was also dis cussed for the neutrino emission and the neutron 
decoupling (Meszaros & Rees 2000b; Paczvhski & Xu 1994^ 
iDerishev et al. 1999; Fuller etal. 2000; Beloborodov 200S^ 
However, the combination of these elements leads to a new 
concept of the VHLF fireball arising from the hot photo- 
sphere, which has not been discussed so far, to the best of 
our knowledge. 

It is useful to refer to an interesting analogy with cosmol- 
ogy. The relativistic protons that are not thermalized after the 
fireball dissipation are similar to the relic particles in the Uni- 
verse, in particular, hot relics such as neutrinos and light dark 
matter, which are relativistic at the freeze out. In this sense, 
a "hot relic fireball" attains a VHLF. The dissipation of the 
GRB fireball is also similar to the reheating of the Universe, 
both of which leads to the entropy production after the birth 
of the fireball. We use "dissipation" for both the entropy pro- 
duction at the collisionless shock and at the thermalization of 
relativistic protons, and "thermalization" for the dissipation 
of the relativistic proton energy into radiation. 

3. FIREBALL DYNAMICS: COLLISIONLESS BULK ACCELERATION 
TO VERY HIGH LORENTZ FACTOR (VHLF) 

In this section, we investigate the fireball dynamics to a 
VHLF via collisionless bulk acceleration in detail, extending 
the idea in the previous section. In ^3.11 we first recall that 
the fireball dissipation is strongly suggested by the observed 
spectral relation. We suggest that the dissipation is caused 

' In this paper, we consider the complete merger case. If we consider the 
reverse and forward shock structure, the further acceleration is seemingly re- 
peatable even after the coasting since the energy can be transfeired from the 
rear to the front shell. However, this is not likely as the fro nt sh ell reacceler- 
ates before completely receiving the rear shell energy [see j|5.4l and Eq. )84H . 



by the mass loading, not by the magnetic reconnectio n or 
neutron decay, for the radiation-dominated fireball. In 33.21 
we use a simple two-body collision to describe the dissipa- 
tion, properly taking into account the relativistic hot motions 
of protons before thermalization to discuss the collisionless 
bulk acceleration. T hen, after giving the photospheric and pi- 
onospheri c rad ii in §3.31 we derive the final coasting Lorentz 
factor in 33.41 In 33.51 we examine thermalization processes 
[p'^ (Bethe-Heitler and photomeson). Coulomb, and plasma 
interactions] other than pp collisi ons, which are relevant in 
some circumstances. We devote 33.61 to consistency checks 
with previous works, and 33.7l to remarks on the connections 
between the collisionless bulk acceleration and the magnetic 
acceleration in making the VHLF fireballs. 

3.1. Fireball dissipation suggested by £pe^]^-L Yonetoku 
relation 

The photosphere model has many advantages for interpret- 
ing the GRB prompt emission (see In this model, we 
identify the spectral peak energy Speak with the fireball pho- 
tospheric temperature T . The fireball is likely radiation- 
dominated since the radiative efficiency is high in most GRBs. 
Under these assumptions, the fireball dissipation is strongly 
suggested by the observed spectral relation, i.e., the e^si^-L 
Yonetoku r elation in Eq. (Ht, as po inted out b y Tho mpson et 
al. (2007) (iThompson et aLll2007t IPe'er et al.ll2007h (see be- 
low). 

In the usual picture that the engine is an accreting black 
hole or possibly a rapidly rotating magnetar, the engine size ro 
is essentially constant about a couple of Schwarzschild radii 
ro ^10^ cm for a black hole of mass Mbh ^ IOM0. However, 
this picture (ro ^ const) leads to a different relation for the 
isotropic luminosity, 

L = AT:rlcaT^ (xT"^, (8) 

from the observed Yonetoku relation L oc T^, where the ob- 
served temperature preserves the initial temperature T ^ To 
for radiation-dominated fireballs. Therefore, the fireball ra- 
dius is most likely reset by the fireball dissipation. Note that 
the relation tracks L oc T after the fireball becomes matter- 
dominant. Although we might be able to transform LocT* to 
L cx by using the matter-dominant track, this is not likely 
since the low-luminosity region becomes radiatively too inef- 
ficient. 

The dissipation takes place at a radius that is much larger 
than the engine radius ro. We may estimate the radius of the 
dissipated fireball at the base of the flow by using the black 
body relation L = 4TT(rph/Tph)^caT* as 

ri = ^ 1 X 10** cm Lgf Teoiev > '"o lO' cm, (9) 

i ph 

where rph is the photospheric radius, Fph is the Lorentz fac- 
tor of the radiating flow, and we have r^ oc L"'/^ if we also 
combine the Yonetoku relation in Eq. ([T]i- 

In general, the actual dissipation radius r,„ ^ T,„ri, is larger 
than the base radius r/,, because, in order to be observed, a fire- 
ball has to have a rel ativistic bulk Lorentz factor T,„ after the 
dissipation (see 33.21 Note that, by definition, the Lorentz fac- 
tor is unity as pulled back to the base of the flow with F cx r). 
The dissipation radius r„, may be comparable to the size of the 
progenitor star r,„ - T„,rh - 10"'-10" cm if r„, - 10^-10^ 
Actually, such dissipation is suggested by numerical simula- 
tions as the jet interacts with the progenitor star (.Lazzati et al] 
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I2OO9I: iZhang etani2004t iMizuta et al.ll200l lIoToh (see also 

§E3]and|5ail. 

In addition, the dissipation has to be associated with the 
deceleration of the fireball. For an impulsive dissipation like 
shocks, the bulk Lorentz factor has to be decelerated by ^ 

— 1 /2 

fb/ro^ 10 L53 times at the dissipation, and for a continuous 
dissipation, the total change of the bulk Lorentz factor would 
be of the same order, because F cx r for radiation-dominated 
fireballs (see Fig. |2]). For the fireball to be decelerated, the 
dissipation seems to be caused by the mass loading, probably 
the baryon loading, not b y the magnetic reconnection or the 
neutron decay (see also § ^5.3l and l5.4l ). 

3.2. Initial condition of dissipated fireballs 

As we have discussed in the previous section, the Epeak-^ 
Yonetoku relation suggests the fireball dissipation under the 
photosphere due to the baryon loading. The main features of 
the dissipation can be described by a simple two-body col- 
lision. In contrast to the previous calculations, we approxi- 
mately take into account the relativistic hot motions of pro- 
tons to discuss the collisionless bulk acceleration. 

We consider a rapid shell that is radiation-dominated with 
internal energy E'^ (oc r"') and Lorentz factor (cx r) ^ 1, 
merging with a slow mass M, with a Lorentz factor Fj. Here, 
Fj = 1 for the matter located near the progenitor star, and F, > 
1 for internal shocks. The total energy in the lab frame is also 
rewritten as 



(10) 



using the duration f,, and the geometry-corrected jet luminos- 
ity Lj = L{9j/2)^. We assume no baryon in the rapid shell 
for simplicity. Although radiation may be already decoupled 
from the baryon-poor rapid shell before the merger, the radia- 
tion is trapped again by the merged shell if the radius is under 
the photosphere. 
The energy and momentum conservation gives 

E'^Tr+M^c^T, = {r,„,Mc^+E;„) F,„, (11) 



£;^F2- 1 +M,cyF2- 1 = (F„„M,c2 + £;) y'F2 - 1(12) 

where E^, and F,„ are the internal energy and the bulk Lorentz 
factor after the merger, respectively, and 



r„„ = F,„F,-./F2-l./F?-l 



(13) 



is the relative Lorentz factor between F,„ and F,. The point 
different from the previous studies is that we take into account 
the prethermalized relativistic motions of protons with F,„j in 
Eqs. (fTTT i and ( fT2] i. Of course, this treatment is approximate 
but valid for order-of-magnitude estimates. If we consider in 
the rest frame of the shocked region, the protons run into the 
shocked region with ^ F,„5 and are isotropized by magnetic 
fields. Here, weak magnetization is sufficient to trap protons, 
and the magnetic fields could be generated by the plasma in- 
stabilities or could be advected from the central engine. Since 
the shock is collisionless for protons before thermalization, 
the proton velocities are not altered so much in the frame of 
the shocked region. Therefore, the random Lorentz factor of 
protons would also be about F,„j in the shocked region. Be- 
fore the complete deceleration, the relative Lorentz factor be- 
tween the preshocked protons and the shocked region is larger 
than F„„. Thus, a fraction of protons would have the random 



Lorentz factor larger than F,„5. However, at least half of the 
protons are shocked after the shocked region is well decel- 
erated, thereby having the random Lorentz factor of ^ F„,j. 
We do not consider the particle acceleration at the shock for 
simplicity. 

We can solve two equations, Eqs. (fTTT i and (fTSI i. for two 
unknowns, F„, and Ei, as 



£';f,+m,c2f. 



(14) 



e;,= J E'; +My + 2E'^,c^T,,-r„„Msc\ (15) 



where F„ = F^F., - yT^-TyFj-T is the relative Lorentz 
factor between rapid and slow shells, and F,, ^ Tr/IT, for 
F,> 1. 

We are now considering the following case: 

(i) the energy is dominated by the rapid shell, E'Tr ^ 
MjC^F,, so that ^'F^ » M^c^ and also the merged shell 
is still relativistic F,„ ^ 1, i.e., F,„j ^ F,„/2Fs, 

(ii) the rapid shell decelerates effectively F„, < Fr/2, con- 
verting the kinetic energy into the internal energy, so 
that 2M,c2f„ > £■;.. 

That is, the slow mass is in the range, 

io-"m0 L5302^„,_3F-2f2 « M,F, « 10-5 Mq L530 'MJe) 

with Eq. ( fTOl l. and thereby, only a small fraction of the pro- 
genitor mass is sufficient to decelerate the fireball efficiently. 
Then, we can simplify Eqs. (fT4] | and (fTsl l as 



1/2 



(17) 



E' 



'^,„,M,c^ (18) 

where we use Eq. ( fTol ) in the last equality in Eq. dTTl ). and in- 
troduce a dimensionless entropy (the radiation-to-baryon ra- 
tio) of the fireball after the merger as 



L 



(19) 



Equations ([TtI i and (fTsT l have important impUcations for the 
initial condition of the dissipated fireballs: 

• First, we can interpret Eq. (fTsT l that the proton energy is 
boosted to a value comparable to the radiation E'^^ and 
almost all the proton energy F^jM^c^ is carried by their 
relativistic hot motions with the Lorentz factor 



7„ 



2F,' 



(20) 



• Second, if the relativistic hot motions are converted into 
the bulk motion, the final bulk Lorentz factor achieves 
the saturation level. 



mip 2 ' 



(21) 



with Eqs. ( [TtI i and ( |20] |. so that the kinetic luminosity 
becomes comparable to the total luminosity, F^Mc^ ^ 
L/2. 
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• Finally, the fireballs should be relativistic even after the 
dissipation in order to produce the observed GRBs; oth- 
erwise, the final Lorentz factor cannot be sufficiently 
high, Fc > 10^-10^, to avoid the compactness problem. 

Note that the bulk Lorentz factor after the merger r„, is 
the same as Eq. (fTTT i even if we omit the new factor r,„j in 
Eqs. (fTTT i and (flSl i as in the conventional case. 

3.3. Photosphere and pionosphere of fireballs 

In this section, we summarize several important radii of 
the fireba lls, in particular, for the photosphere and pio- 
nosphere (iRees & MeszarosllIOOl: lM6szaros & Reesll2000i: 
iMeszaros et al.ll2002l) to derive the final coasting Lorentz fac- 
tor in the next section. We assume that the fireballs are created 
intermittently with a variability timescale f,, ~ r/,/c, which is 
determined by the size at the base of the flow, for simplicity, 
although we may apply the following arguments as long as 
the timescale is longer than the causal timescale, f,. > r/cT^. 

Baryonic photosphere rph is determined by rj- = 1, where 
the Thompson optical depth to electrons associated with pro- 
tons, Tj = n'^ajr/T, is unityQ The Lorentz factor is F = 
[r/rh^r]}, the comoving width of the shell is A' = [r, r/,?7], 
and the comoving volume of the shell is V' = Aur^hJ = 
[Anr^ ,4TTrirbr^] when the radius r is [< rbri,> ri,rf\. Then, 
with the baryon density n'^ = Lty/r]mpC^V' , we can derive the 
photospheric radius as 



_{TfJ^'q^l^ for rj > T]^ 
fb Xvtv'^ for V<V*, 



rph 



where 



\AiimpC^rb 



1/4 



1 X lO^L^/V,,, 



1/4 



(22) 



(23) 



is a critical dimensionless entropy. We can apply the above 
relations to the wind regime, r > ct^.rf' ~ r^jTy^, where the suc- 
cessive shells expand their thickness and overlap through in- 
ternal shocks, because n'^ = L/A-Trr^mpC^rjT and the Lorentz 
factor is saturated at r = r;,ry with T = 1]. 

Pionosphere is defined by Tpp = 1, where the optical depth 
to pp collisions (i.e., pionic optical depth) is unity. Repeating 
the previous calculations, we have 



rb 



where 



'npp- 



( 



4 -3 



1/4 



for 77 > 'qpp 
for r]<ripp, 



(24) 



\4T:mpC^rh 



^pp 



1/4 



'50QLfrlfi25) 



An photosph ere can be formed beyond the bary- 
onic photosphere (Pilla&Loeb' '1998'; 'Guetta et al.' '20011; 
iKobavashi et al. 2002; Meszaros et al. 2002; Li et al. 2003; 
iRees & Meszaros 2004; loka e t al.ll2007} . Although there are 
some uncertainties in the amount of e^, the actual abun- 
dance will be between the following three cases: 

(1) The e^-p equal mass case, in which the rest mass energy 
density of is equal to that of baryon. We can obtain the 



2 In Ref. IMeszaros etaTI UM^ . TT = n' ajh! was used in the discrete 
shell regime. However, we think that tj — n' ajrjV is correct since photons 
can travel only a distance ~ r/F within the comoving time at radius r. 



photospheric radius r± by replacing nip and nig in the baryonic 
photosphere case as 



[^_\r\j^r\ 'Z-' for rj > ri± 
n, \r]±'n~^ for V<V±, 



(26) 



where 



Lot 



)'''=(^)"'.7*-7xlO^^L;/Vf27) 

(2) The feasible case, in which the density is determined 
by the balance between annihilation and creation, where 
e^ is created by pp collisions. This is guaranteed in the pres- 
ence of relativistic hot protons, which remains not thermal- 
ized after the fireball dissipation. 

The inelastic cross section for pp collisions is about app ^ 

3 X 10"^* cm"^ above the pion production threshold ^ 140 
MeV. The inelasticity Kpp is ^ 0.5, so only a few collisions are 
required to extract most of the energy of the primary particles. 
The energy is initially given to produce tt* and tt", where the 
TT multiplicity is typically ~ 1-3 near the threshold ^/s ^ 
1 GeV and w eakly depends on the center-of-mass energy as 
Mtt cx In-y/s (iGrosse-Oetringhaus & Revgersll2009D . The pi- 
ons immediately decay as tt""" — ^ 11^ + — >■ e'^ + Ve + Vp. + i^^, 
7r~ — > /i~ + — > e~ + + z^p + P^, and vr" — > 7 + 7, and the 
decay gamma-rays create e^ via 77 interactions. Thus, the 
minimum multiplicity without considering the following 
cascade is 



1/4 



AmneC rt, 



2M, 



(28) 



Alternatively, since the injected has a large Lorentz fac- 
tor 7j_ ~ Kpp^'pmp/2M..^nig, the electromagnetic cascade will 

follow and could achieve the maximum e^ multiplicity. 



' Kpp-fpmp/nie, 



(29) 



where 7^ is the comoving Lorentz factor of protons. The min- 
imum e^ multiplicity would be appropriate when the syn- 
chrotron cooling dominates, since the synchrotron photons 
are typically soft < 1 MeV, while the maximum multiplic- 
ity would be valid when the inverse Compton cooling dom- 
inates, since the relativistic e^ can scatter photons to high- 
energy ^ 1 MeV. Hereafter, we parametrize the multiplic- 
ity as 

M±=f±j'pmp/me. (30) 
Then, equating the annihilation rate 



3 

nj_ = —n^n_<jjc 



(31) 



with the creation rate 

n'^ = M±n'^pappC, (32) 

we obtain the feasible density n± = n+ ^ n- as 

1/2 

, -f±l'p 1 

3 ar me ' 



(33) 



which can be much larger than unity. 

(3) The maximum case, in which the comoving radiation 
energy is almost converted to the rest mass energy. This 
extreme case might be reahzed if the radiation spectrum has 
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Fig. 3. — r]-r diagram to re ad normalized radii r/ri, an d th e final coasting 
Lorentz factor Fc in Eq. (39) (red thick line) and in Eq. (47) (with efficient 
Bethe-Heitler processes; red thick dotted line) as a function of the dimension- 
less entropy 'q = LlM(r (the radiation to baryon ratio) . A fireball can achieve 
a VHLF, Ff > >7, , if the dissipation radius r,,, in Eq. (53) exceeds the piono- 
sphere r,,,, in Eq. (22) (pionic optical depth Tpp = 1). The maximum coasting 
Lorentz factor reaches Fc.niax ^ 10* in Eqs. (41) and (42). We also show the 
baryonic photosphere Tph in Eq. (24) (Thompson optical depth tt = 1), the 
photosphere r± in Eq. (26) for the e^-p equal mass case, and the satu- 
ration radiusrj/£^= T). C ritic al entropies ?7«, 77^1, 7742, and t])^ are given by 
Eqs. (23). (36), (^, and (38), respectively. 

a significant fraction of energy above the pair production 
threshold and the pairs develop an electromagnetic cas- 
cade. With the density = LtyjVmeC^V' , the optical 
depth condition tj = n'^ ^^-^<7Tr /T = 1 gives 



'"±,1 



rb 



for 
for 



V>V± 

V < v±- 



(34) 



The rest mass energy surpasses the baryonic one, be- 
cause the proton energy, comparable to the radiation energy 
in Eq. (fTsT l. is carried by relativistic motions, not by the rest 
mass energy. The acceleration terminates at this radius r± 
since the rest mass energy almost equals the total energy. 

3.4. Final coasting Lorentz factor of dissipated fireballs 

Figure [3] shows the 77-r diagram, dRees & Meszaroslll994t 
iMeszaros & Reesll2000albii) which is useful to read out the fi- 
nal coasting Lorentz factor. In the conventional picture, we 
usually start the fireball evolution from the base of the flow, 
r/rb = 1, for a given entropy 77. The fireball expands with 
r cx r as far as below the saturation radius and the pho- 
tospheric radius rph. Then, the maximum coasting Lorentz 
factor is ^ 77* in Eq. ( l23T l. 

However, the starting radius is different from the base, 
r/rb = 1, for the dissipated fireball case because the observed 
fire balls have to be relativistic after dissipation, as discussed 
in ^3.21 Then, using F oc r, we can plot the starting radius r^ 
in the diagram (dotted line in Fig.|3]l as 

(35) 
rb 

where F,„ is the bulk Lorentz factor after dissipation in 
Eq. ([TtI i. derived from the energy and momentum conserva- 
tion in Eqs. (fTTT i and (fT2] ). and Fig. [3] shows the case F, = I. 
Here, we take the independent model parameters as the di- 
mensionless entropy 77 in Eq. fT% . the Lorentz factor of the 
slow mass F^, the total isotropic luminosity L, and the engine 



size ro (< r;,) in ^3.11 We can derive the other parameters, 
such as the observed temperature T in Eq. ([T]i, the radius at 
the base of the flow ri, in Eq. (|9]l, the dissipation radius r„, 
in Eq. ( |35] |, and the Lorentz factor of the rapid shell before 
merger F, = r,„/ro. 

As long as the dissipation radius is below the pionosphere, 
rm < rpp, i.e., in the low to moderate entropy range 77 < rj^i 
(see Fig. O where 



Vk. = rfHi' 



2/5 



(36) 



with Eqs. (l24l i and (l35T l. the fireball evolution is essentially 
similar to the conventional case. This is the case (A) in 
Fig. [3] The relativistic hot motions of protons, acquired at 
the merger with a radiation-dominated fireball, are quickly 
thermalized into radiation via pp collisions under the pio- 
nosphere. The dissipated fireball goes back to the standard 
radiation-dominated fireball, and hence, the maximum coast- 
ing Lorentz factor is again ^ 77, in Eq. ( l23l) . as the radiation 
escapes from the photosphere r = rph. Such a fireball is some- 
times called a reborn fireball (Ghisellini et al. 2007). 

The evolution is completely different if the dissipation ra- 
dius exceeds the pionosphere (see ^3.5l for other thermaliza- 
tion processes). This is the case (B) in Fig. [3] Relativistic 
hot motions of protons are not effectively thermalized, which 
are converted into the bulk motion, leading to a VHLF up to 
the saturation value F^ ^ 77 in Eq. ( I2TI ). Here, the acceler- 
ation r c>c r continues even beyond the photosphere because 
the pressure is provided by the collisionless hot motions of 
protons. For a large dissipation radius exceeding the bary- 
onic photosphere r,„ > rph, the creation is necessary to trap 
the radiation that boosts the proton ener gy. A lthough the 
abundance is somewhat uncertain (see 33.31 ), we first derive 
the upper limits on the coasting Lorentz factor by considering 
the most favorable case. For a dissipation radius r^ > r± in 
Eq. ( |26] |, i.e., in the range 77 > rik2 where 



m2 = r-y'vT 



1x10" FT^M^ 



-2/5 
b.S ' 



(37) 



the e^ rest mass necessary for trapping the radiation domi- 
nates the baryon rest mass, suppressing the coasting Lorentz 
factor below the saturation value rj in Fig. [3] An upper limit 
to the dissipation radius is r,„ < r±.max in Eq. ( l34l ). i.e., in the 
range 77 < 77^3 where 



:4xl0^r;'L^fr-f, 



(38) 



because the fireball larger than r± is too rarefied to trap the 
radiation even if e^ is maximally created. Thus, a dissipated 
fireball with 77 > 77^:3 is never formed via physical processes. 
Since the rest mass equals the total energy at r„, = r± max 
in Eq. ( |34| ). the dissipated fireball terminates its acceleration. 

The final coasting Lorentz factor (red thick line in Fig. O 
may be summarized as 



for 77 < 77* 

for 1]^ <ri < rjki 

) V for r]n <'n < Vki 

I Vk2('n/ rik2)^ for r]k2<ri< Vk3, 



(39) 



where the index A is determined by the condition Tc = r„,/rh = 
V^sm at 77 = 770 (see Fig.O as 



A = - 



, / j^-l/2 -1/2 

ln(^77i2F, 77^3' 
In {Vk3Vk2) 



(40) 
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TABLE 1 

Final shares of the initial total luminosity L as a function of the 
radiation-to-baryon ratio 7y. the lumi nosity is shared by the photospheric 

RADIATION Lph, NEUTRINOS (SEE ji5.lt . AND PROTONS AND IN THE KINETIC 
FOR M =^P+L±' WHICH IS LATER DISSIPATED VIA INTERNAL/EXTERNAL SHOCKS 
(SEE j34.U . In 34.11 WE CONNECT Lph WITH THE BAND SPECTRAL COMPONENT AND 

WITH THE EXTRA HIGH-ENERGY POWER-LAW COMPONENT (PL IN SHORT), 
RESPECTIVELY. NOTE THAT THE COASTING LORENTZ FACTOR Fc ACHIEVES A VHLF IN 
THE HIGH ENTROPY RANGE r]iii < rj [SEE FIG.[3] aND EQ. ll39tl. CRITICAL ENTROPIES 
V*' Vkl, Vkl, AND ria ARE GIVEN BY EQS. JIbQ?^, ( [371 , AND ( [38t . RESPECTIVELY. 



ri Tc Lph [~Band] L/, [~PL] Spectrum L,, 



1 <ri Krjt ^ 
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V* <V< Vkl ' 
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<L 


Band 


~L 


Vkl <V< Vkl • 


10' 


V 


~L 




Band+PL 


<L 


Vk2<V< Vk3 ■ 


- 107 


Vkliv/Vkl)^ 


~L 




Band-HPL 


<L 



The appropriate coasting Lorentz factor for 77,^ < 77 < 77*1 

would be 77, instead of the slightly smaller rfj^rf^l^ in 
Eq. ( [39] l because most of the electrons above the photosphere 
can still scatter with a decreasing fraction of free-streaming 
photons, and continue accelerating as long as the comoving 
Compton drag time t'^^^^ = nipC^ /ca rU^ is less than the co- 
moving dynamical time fj^ = r/cF (iGrimsrud & WassermanI 
ll998tlMeszaros & Reesll2000ah . The ratio of these two times, 
4agAdyn = '^T^rripC^ rT^ / Lar = (r/ ij^rb)'*, gives the coasting ra- 
dius where the radiative acceleration ceases at r/rt = rj^, > 
rph/n for r] > rj^,. 
The maximum value of the coasting Lorentz factor is 

r.,™ax = %2 - 1 X 10^ T:y'L%'r;f, (41) 

(see Fig. O, which is a VHLF, much larger than the conven- 
tional upper limit 7;, ^ 10^ in Eqs. (|5]l and ( l23T l. This max- 
imum value is realized by the e^-p equal mass case when 
the rest mass energy is equal to the baryonic one. If the 
actual abundance is less than this case, the radiation es- 
capes without completely transferring its energy to the proton 
component, leading to a smaller Fc niax- A more conservative 
estimate of the maximum Lorentz factor is provided by the 
feasible case of the creation in Eq. ( [33T l of ^3.31 which 
gives the feasible optical depth to trap the radiation at the dis- 
sipation. By solving tt = 2n'j_f77-r„,/Fm = 1 in terms of rj, 
with Eqs. (EU, d^, and nj, = L/A-TrrlmpC^rjT^, 

I 



we may identify its solution with the most conservative max- 
imum Lorentz factor as 

pc _ 16 (Tpp nip p-7/9^16/9 

^^''^--[ 3 ar m/^J ' ^* 

^SxloYfr-y^LfrlT, (42) 

which is still a VHLF unless f± is very small in Eq. ( |30l l. 
Therefore, we conclude that the coasting Lorentz factor can 
attain a VHLF, F^ ^ 10^-10^, if the dissipation radius is above 
the pionosphere, i.e., in the high entropy range with rj> rjk\^ 

lO^Fj^^^ (see ^3.51 for other thermalization processes). It is 
interesting to note a sharp rise of F^ at 77 = 77^1, where a slight 
change of the baryon loads could lead to a large difference in 
the coasting Lorentz factor AFc/F^ > 10. 

The initial fireball energy is shared by the photospher ic ra - 
diation luminosity Lph, the neutrino luminosity Li, (see 35.11 ). 
and the kinetic luminosity of protons and e^, = Lp + L±, 
where the kinetic lumino sity is later dissipated via internal 
and external shocks (see ^4.1t . Table [T] summarizes the final 
shares of the initial total luminosity L as a function of 77. As 
we discuss in 34.11 if we connect Lph with the Band spectral 
component and with the extra high-energy power-law com- 
ponent, the equal contributions from both components could 
suggest the VHLF range with 7;ai < r] in Table [1] Fig. |3] and 
Eq. (|39] l, and the GeV onset delay phase could suggest the 
moderate entropy range 7/, < 77 < 77^:1 (see 34.4l l. 



3.5. Thermalization: pp, pj. Coulomb, and plasma 
processes 

We have considered pp collisions for the thermalization 
process of relativistic protons that are entrained in the fire- 
ball dissipation. In this section, we also examine other pro- 
cesses, pj (Bethe-Heitler and photomeson) and Coulomb in- 
teractions, which are relevant in some circumstances. (The 
bremsstrahlung emission is not effective for fiducial parame- 
ters.) 

The main competing interaction is Bethe-Heitler photopair 
production (pj pe^e~). The threshold energy of photons 



normalized by the peak energy is about 

ImeC^I-l'p 



-peak 



-2F, 



3 L-y^F, 



,(43) 



where r > r„, and we consider the adiabatically cooling pro- 
tons with the comoving Lorentz factor 7^ ^ (r,„/2Fj)(r,„/r) 
in Eq. (|20l), photons with energy e^^^j. = (r/,/r)£peak in Eq. ([B, 
and the dissipation radius r,„/r/, ~ F,„ in Eq. ( [35] ). The ratio 
in Eq. ( 1431 ) is usually larger than unity. Therefore, the Bethe- 
Heitler process is not effective for a thermal photon spectrum 
with a cutoff above the peak energy Epeak- 
Whereas, if the photon spectrum is nonthermal, the protons 
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can interact with high-energy photons to produce pairs. For a 
typical Band spectrum, the protons can cool via Bethe-Heitler 
processes even above the pionosphere without pp collisions. 
The Lorentz factor of protons decreases to a value 7'™ for 
which the Bethe-Heitler cooling optical depth is about unity, 



1, 



(44) 



where we estimate at the dissipation radius r = in Eq. (l35T l. 
and approximate the photon number density above the thresh- 
old as 



«;(<>2m,cV7T) 



(45) 



with the high-energy photon index /3 ^ 2.5, the cross section 
(Jbh ^ (e^/^c)(3 /87r)CTr and the inelasticity Kbh ^ me/mp for 
the Bethe-Heitler process. We note that the value KbhCbu ^ 
5 X 10"^' cm^ changes by less than a factor of 3.5 for the 
range 5 < v'^ip/nieC^ < 10^ jChodorowski et al.ll 19921) . From 
Eqs. (l44l) and ( |45] ), we obtain the Lorentz factor of protons 
after the Bethe-Heitler cooling as 



7, 



,BH 2meC^T„ 



-peak 



p3/2^3/2^-4 Speak (Tt(J3-1) 

nipC^ ^^bho-bh 



(46) 



with Eqs. (|23]) and ([35]l and Tms'^.^^- = Cpeak, where 7'™ 
does not drop to unity in most cases, differently from the 
pp collisional case. We can neglect the Bethe-Heitler cool- 
ing above the dissipation radius as the number density of 
photons above the threshold rapidly decreases as n'^{v'^ > 
2m(.c^/7'™) cx r"^2;3+i) leftover relativistic motions of 

protons in Eq. ( |46l ) are converted into the bulk motion via 
coUisionless bulk acceleration. Therefore, in the case of the 
effective Bethe-Heitler cooling, we derive the final coasting 
Lorentz factor as 



^ c 



,BH 



where 



-) 

2(/3-l)/3 



3/2(/3-l) 



mpC^ ^^bhO-bh 



(47) 



-,2/3 



Epeak (Jt(J3- 1) 



^6x 10" r 



-(2/J+I)/3,2/3 -2/3 2(/3-2)/3 
; -^53 'fc,8 ^peak.MeV 



(48) 



with Eq. (I23]l. In Fig. |3] we plot L^" of Eqs. (|47]i and 
(l48T l. We can see that the Bethe-Heitler process can reduce 
the proton kinetic energy by r^/Lc r^"/r/ 0.03- 
1 times (at most) in the range of the dimensionless entropy 
Vki < 1] < TjBH- Nevertheless, the coasting Lorentz factor can 
attain a VHLF even if the Bethe-Heitler cooling is most ef- 
fective. The actual evolution would be between and 
in Fig. |3] but is difficult to evaluate exactly since it depends 
on the photon spectrum below the photosphere, which could 
be thermal or could already be the Band spectrum (see 35.2| i. 
We also note that the coasting Lorentz factor could be larger 
than if the pairs produced by the Bethe-Heitler pro- 
cess make a photosphere above the coasting radius since the 
radiative acceleration continues up to the photospheric radius. 

The other photoprocess is the photomeson interaction, 
which produces one or more mesons, mostly pions. The 



threshold is higher by niTr/nie ~ 280, while the mean cross 
section Kpm<Jpm ^1 x 10"^ ^ cm^ is ^ 140 tirnes hig her than 
the Bethe-Heitler process (IChodorowski et al]ll992l) . Then, 
the ratio of the photomeson cooling to the Bethe-Heitler cool- 
ing is 



^BHCBH 



1-/3 



'3x10-2 for (3 = 2.5, (49) 



and ^ 0.5 for (3 = 2. Thus, the photomeson cooling is sub- 
dominant for a typical Band spectrum with (3 > 2. 

The Coulomb collisions with also contribute to the pro- 
ton cooling because the temperature of is usually kept at 
a nonrelativistic value, m uch below the proton temperature, 
by the Compton cooling (IBeloborodovll2009 ) . The effective 
radius of the Coulomb cooling is determined by 



n±- — ajlnA^ 1, 
1 nip 



(50) 



where InA 10 is the Coulomb logarithm. That is, the 
Coulomb cooling is effective if the Thompson optical depth is 
larger than tj > nip/nielnA ^ 200. Assuming the feasible 
density in Eq. ( |33] ) at the dissipation radius r,„ in Eq. (|35] |. we 
can estimate the range of the dimensionless entropy rj where 
the Coulomb collisions are effective as 



V<Vc'- 



: „16/V-7/9 
' '/* ^ s 



3 ar mp 



2/9 



4 7-4/9 -i/9p-7/9. 2/9 



^6xio^l;^%, 



± 



(51) 



with Eq. d20] i and the relation «Jj = L/Airrl^mpC^riTyn. Then, 
for the maximum multiplicity f± ^ 0.5 in Eq. (|30] |. we 
have rjc > rjk\ ^ 2x lO'* in Eq. ( |36] | and Fig. [3] that is, the 
protons cool via Coulomb collisions for rjki < rj < rjc even 
above the pionosphere without pp collisions. (Note that the 
pp collisions dominate the Coulomb collisions for 77 < ?7jti 
even in this case). The random Lorentz factor of the protons 
decreases, reducing the coUisionless pressure, although it may 
not completely vanish since the density also decreases as 
the protons, i.e., the energy source of e*, cool down. Then, 
the coasting Lorentz factor for 77^1 < 77 < 77c could not at- 
tain the saturation value Fc = in Eq. ( |39] | and Fig. [3] Never- 
theless, the ratio ijc/Vki ^ 3 is just a factor and does not de- 
pend on the parameters so much {rjc/rjky oc L^/'^^r^^'^^^Fj'^^). 
Thus, the critical entropy 77^1 for the pp collisions in Eq. ( [36] l 
is still a good indicator of the VHLF fireball formation. 

In addition to the above two-body processes, the collec- 
tive plasma processes could dissipate the relativistic proton 
energy. At least, a fraction (eg ^ 0.01) of the proton en- 
ergy could be converted into the magnetic field via the Weibel 
instability (Medvedev & Loeb 199 9D or the macroscopic tur- 
bulence ( Goodman & Ma cFadyenI l2008t ISironi & GoodmanI 
i2007; Zhang et al., 2009) . However, the conversion into the 
magnetic field does not reduce the coasting Lorentz factor be- 
cause the magnetic field also provides the relativistic pressure 
to expand the fireball (see ^3.7b . Rather, it opens up a way 
to the VHLF since the magnetic pressure can persist even af- 
ter the proton cooling via two-body processes. On the other 
hand, the dissipation of the proton energy into electrons is 
harmful for making the VHLF. However, the ene rgy conver- 
sion fraction is usually less than half (e^ < 0.5) dSpitkovskvl 
I2OO8I: iToma et aTl 120081) . The plasma dissipation would be 
also suppressed after the protons are isotropized. Therefore, 
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the plasma processes are unlikely obstacles to the VHLF for- 
mation. 

3.6. Consistency with previous works 

We have discussed a physical mechanism of the collision- 
less bulk acceleration to create VHLF fireballs with F^. > ry* ^ 
lO-' up to Tc < 10* in the previous sections. Since such a 
VHLF is somewhat extreme and has not been considered se- 
riously, we examine whether a VHLF is allowed by the previ- 
ous observations. 

(1) A high-energy cutoff due to creation pro- 
vides in formation on the bulk Lorentz factor of a 
fireball (iBa-ing & Hardi ngI [T997t Oth wick & Saril l200ll 
iRazzaque et 30 12004; M urase & Iok"al l2008; Aoi et al. 2009). 
This method usually gives a lower limit on the bulk Lorentz 
factor, which is about Tc > 10^-10^ so far (Abdoetal. 



|2009ba, 2010c) and consistent with a VHLF (see S43\fov 
GRB 090926). Further observations are needed, and hope- 
fully, more elaborate observations are necessary because the 
exponential cutoff is usually smoothed to a broken power- 
law by multizon e effects, wh ich make it difficult to iden- 
tify the cuto ff (Aoi e t al.||2009t [Lil l201(ji: iGranot et al.ll2008t 
iBosniak et al. 2009) . 

(2) For the internal shocks to take place before the external 
shock, the minimum Lorentz factor in the flow has to be below 
r < 3 X 10". (Rees & Meszaros 1994; Sari & Piran 1997a) 
This is because the internal shock radius becomes larger for 
higher F, 



rsh = 2r^ct, - 6 X 10'^ cm F?f„. 



■3- 



(52) 



On the other hand, the afterglow starts when the reverse shock 
crosses the shell. Since a VHLF evolution is in the so-called 
thick shell case, the reverse shock crosses at the duration time 
T. After that, the hydrodynamic evolution enters into the self- 
similar phase, which may be described by the adiabatic condi- 
tion E ^ {Ait /Tjj^r^nmpC^ . Eliminating 7 with T ^ r^^/AT^c, 
we obtain the reverse shock crossing radius as 



■7x WcmE. 



1/4 
53 



-1/47,1/4 



(53) 



which do not depend on F. We note that a VHLF with 
F > 3 X 10"* is allowed if it is decelerated by internal shocks 
with slower shells before external sh ocks. This is even fa- 
vorable for effic ient inter nal shock s (Koba vashi et al.l [19971 ; 
iKobavashi & Sa ri 2 001; Beloborod ov 2000), and also for the 
GeV onset delay in ^441 

(3) By identifying the peak time of the afterglow light curve 
with the decelerati ng time of the eject a , we can con strain 
the Lorentz factor dSyr& PiranI 1 1999bt iRvkoffet all 120091: 
iGhirlanda et ani20ia iLiang et alJl2009 ). Several results do 
not imply VHLF values but typically F ~ 100-600. However, 
this method is only applicable to the so-called thin shell case, 
whereas a VHLF evolution is likely a thick shell case. In the 
thin shell case, the peak time of the afterglow light curve is 
expected to be later than the prompt emission, while in the 
thick shell case, the peak time is comparable to the prompt 
duration and difficult to observe. This method also gives a 
lower limit on the Lorentz factor if the flow is decelerated by 
internal shocks with slower shells before external shocks. 

(4) Reverse and forward shock emission in the early after- 
glows constrai ns the Lore ntz facto r via spectral and tempo- 
ral modelings dSari & Pir an 1999a; IZhang et al.ll2003l) . Sev- 
eral results do not imply the VHLF. However, the early after- 
glow modehngs are confronted with difficulties to interpret 



t he s teep a nd shall ow decay discovered by Swift (see also 
EH (Zha ng etalJj 2006; lok a et alJ |2006; Panaitescu et aH 
I2006t .Zhang et alJ|jBo6^,HuangetaLll2006i;.Sato et al.i,200C 
This method also gives a lower limit on the Lorentz factor if 
the flow is decelerated by internal shocks with slower shells 
before external shocks. 

(5) We can infer the Lorentz factor by identifying the spec- 
tral peak L{v = gpeak) with the ther mal emission component 
dPe'er et al.1 120071 iRvde et alJl20Tol) . Pe'er et al. and Ryde 
et al. applied this method to several bursts including GRB 
090902B and suggested F - 300-800 below a VHLF How- 
ever, this method potentially has two solutions, one of which 
may provide a low Lorentz factor, whereas the other could be 
the real case with a VHLF. To be more precise, this method 
initially discriminates two possible cases, rph > and rph < 
rs, where it is possible to determine the Lorentz factor only 
in the case of rph > rj. This case also provides a consistency 
inequality (rj < rj^,) between observables, so one might think 
that this is the solution if the inequality is satisfied. However, 
it is logically not a sufficient condition but just a necessary 
condition to satisfy the inequality, so that we cannot exclude 
the other case rph < r,, which allows a VHLF. 

(6) Zhang and Pe'er (Zhang & Pe'er 200p; |Fa3 IMol) 
claim that the predicted thermal component is not consis- 
tent with the observation, suggesting that the outflow is not 
radiation-dominated but Poynting-dominated. They estimate 
the maximum temperature allowed by the observation of GRB 
0809 16C as T°l^^^ = 50 keV, using the relation rh = c5t°^ 

with the observed variability time St°^ = 0.5 s, and conclude 
that this is below and contradicts with the observed peak en- 
ergy. However, the central engine size may be smaller than 
cdt°^, so that we can raise the maximum temperature with 
Tph.max oc r"'/^ to fit the observation. In addition, the pho- 
tosphere model usually assumes that the black body spec- 
trum is modified to a Band spectrum by Comptonization via 
magnetic reconnection, neutron collisions or repeated shocks, 
([Thompson et al. 2007; Meszaros & Rees 2000a; lokaetaQ 
12007; Beloborodov 2009; Giannios & S pruit 2007; Gianniol 
12008!) although t he a c tual mechanism has not yet been re- 
vealed (see also ^5.2\ dioka et alJl2007h . Then, it is possible 
to reproduce the observed spectrum by the photospheric emis- 
sion. 

(7) Zou and Piran ( Zou & Piran| |201(il) gave an upper limit 
on the Lorentz factor by requiring that the observed deep min- 
imum in the prompt phase should be above the early external 
shock emission. The results are F < 10^, not a VHLF. Since 
they present an analysis only for the thin shell case, it is de- 
sirable to calculate the thick shell case, which is the likely 
case for the VHLF evolution. This method also gives a lower 
limit on the Lorentz factor if the flow is decelerated by inter- 
nal shocks with slower shells before external shocks. 

In conclusion, the VHLF fireballs are currently consistent 
with previous observations. 

3.7. Connection with magnetic acceleration 

A VHLF can also be achieved by the magnetic accelera- 
tion of a fireball. If the energy density is dominated by the 
magnetic fields, the magnetic pressure expands the fireball 
up to the equipartition of the comoving energy density be- 
tween the magnetic fields and matter (protons and/or elec- 
trons) (Mes zaros & ReeslfT997h . 

The difference between the magnetic acceleration and the 
coUisionless bulk acceleration is in the radiation fraction 
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L^h/L after the evolution. In the collisionless bulk accelera- 
tion, the radiation fraction Lpi^/L is about unity for the VHLF 
fireballs in Table[T] Whereas, in the magnetic acceleration, the 
radiation fraction Lph/i depends on the initial condition (i.e., 
the radiation and fraction) and the magnetic reconnection 
during the expansion, thereby not determined from the first 
principle. 

We can also consider the magnetic component in the colli- 
sionless bulk acceleration. Firstly, the magnetic field could be 
advected from the central engine. In this case, the final radia- 
tion fraction L^h/L becomes less than unity because the mag- 
netic component occupies a fraction of the final luminosity 
as a nonradiative component. The magnetic component from 
the central engine could be relevant for the short GRBs as 
discussed in ^5.71 Secondly, the magnetic field could be gen- 
erated by the relativistic proton com pone nt via plasma pro- 
cesses at the fireball dissipation (see ^3. 5b . This does not re- 
duce the final radiation fraction Lph/L. The magnetic pressure 
could persist even after the proton cooling, expanding a way 
to the VHLF. 

4. FIREBALL SPECTRUM: HOT 
PHOTOSPHERE-INTERNAL-EXTERNAL SHOCK 
SYNCHROTRON MODEL 

In this section, we apply VHLF fireballs to the GRB emis- 
sion. In 34.11 we first discuss the GRB emission site from 
the kinematical point of view without detailed spectral mod- 
elings. In 34.21 we calculate the internal shock synchrotron 
spectrum and find that the VHLF models have advantages on 
providing a single rising power-law spectrum over >7 energy 
digits and also the high internal shock efficiency to supply suf- 
ficient energy to the power-l aw co mponent that is comparable 
to the Band component. In 34.31 we discuss the possible ori- 
gins of the spectral break at approximately 1 .4 GeV observed 
in the extra component of GRB 090926, which was suggested 
as the creation cutoff for F ^ 600. In the VHLF mod- 
els, the spectral break could be the synchrotron cooling break 
for F ~ lO'*, or the maximum synchrotron cutoff, particularly 
limited by the dynamical time, for F ^ 10^. In 34.41 we sug- 
gest that the GeV onset delay reflects the timescale for the 
baryon loading rate to change at the dissipation radius, i.e., 
the environmental change around the progenitor star controls 
the GeV onset delay. We can predict the delay time and also 
its parameter dependences, which are found to be consistent 
with the observations. 

4. 1 . Emission sites 

The GRB emission has at least two, and probably three 
components: (1) the usual Band component (a broken power- 
law with a peak at Cpeak ^0.1-1 MeV), (2a) high-energy emis- 
sion (> 100 M eV) that lasts longer than the MeV emission, 
dHurlev et al.| [T994: Gonzalez et al. 2003) sh owing a powe r- 
law deca y, as revealed by Fermi/LAT ( Abdo et al. 2009c aQ, 
I2010blcl) . (2b) high -energy variable emission (> 100 MeV) 
that looks correlated with the MeV emission with a short rise 
and decay time scale dAbdo et al.ll2009a) . Most of the spec- 
tra are well fitted by the Band function even up to ~ 10 GeV, 
while an additional distinct component at > 10 MeV is also 
found, at least in GRB 090510 and 090902B out of - 12 
Fermi/LAT bursts. This additional component is fitted by a 
single power-law that slightly rises in vFi, cx and of- 

ten extends below < 20 keV over ~ 7 energy digits. (How- 
ever, note that no other experiments have confirmed the low- 
energy extension of the power-law component.) 



Figure |4] shows our model in a schematic way. This is es- 
sentially the same as the photosphere-in ternal shock model 
proposed by Toma et al. ( iToma et ani201 0) (see 34.2l for dif- 
ferences in the Lorentz factor and hence the emission mecha- 
nism). In this model, we consider the following emission site 
for each emission component, 

(1) Band component: Photospheric emission, 

(2a) Long-lived power-law component: External shock, 

(2b) Variable power-law component: Internal shock. 

These assignments seem reasonable from the kinematical 
point of view, even without detailed spectral modelings, be- 
cause of the following reasons: 

• A photospheric origin of the Band component easily ex- 
plains the high radiative efficiency, in addition to the 
low-energy spectral index and the stability of the spec- 
tral £peak relation (see ® . The radiative efficiency (i.e., 
the ratio of the photospheric luminosity Lph to the total 
luminosity L) can be more than ^ 50% if the radiation 
to baryon ratio is in the moderate to high entropy range 
rj> rjf, ^ IQr', as summarized in Table[T]with Eq. (|23] |. 
The dissipation under the photosphere is suggested by 
the Epeak-i- Yonctoku relation (see 33.11 1. 

• An external shock origin of the long-lived power- 
law component has an advantage for explaining its 
longevity and power-law decay. The high-energy 
long-lived emission could be consistent with the sim- 
ple synchrotron emission from the external forward 
shock, wi th either adiab atic (Kumar & Barniol Durar 

Barniol Duran & Kumar 
baoetal. 2009; 



2009a 


b 


2009a 


b 



2010; Corsietal 
De Pasquale et al. 201 
Pandev et al. 2010) or radiative shock ( Ghirlanda et al.l 
2010; Ghiselli ni et al. 2010) and probably flie Klein- 
Nishina effect ( Wang et al. 20Tot iNakar et al.1 E^009). 
Although the maximum synchrotron frequency has not 
been detected (Li 2010b; Piran & Nakar 2010), the 
long-lived component co uld be p roduced by t he in- 
verse Compton emission dZou et al. 2009; Wang et al 
20101: ICorsi et a l. 2009a; Neamus 2010; M urase et al 



2010I) . As we can see from Table [T] the kinetic energy 



of the external shock is comparable to the total energy if 

r^>m^ lO^r;^'''* in the VHLF range and ?7 < 77, ^ 10 
in the low entropy range with Eqs. (l23l l and (|36] |. 

An internal shock origin is naturally invoked for 
the variable high-energy component because the ex- 
ternal shock cannot usually produce fast variability 
([Sari & Piran 1997b; loka et al. 2005). The external 
shocks cannot have the fast decay of the light curve 
(floka et al.. 200 5). although the fast rise could be pro- 
duced by the finite acceleratio n time of particles that ra- 
diate at the observed energy (iBarniol Duran & Kumad 
I2OIO). It is favorable to have VHLF fireballs in addition 
to slow ones, i.e., a large dispersion in the Lorentz fac- 
tor, to raise the internal shock efficiency for converting 
the kinetic energy into radiatio n ( Kobavashi et al. 199'3 
iKobavashi & Saril 120011: iBefoborodov 200 0), because 
the variable power- law component has a comparable lu- 
minosity to the Band component in several bursts. It is 
interesting to note a sharp rise of the coasting Lorentz 
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Fig. 4. — Photosphere-internal-external shock model jToma et ani2010f) 
in which the photospheric emission produces the Band spectrum, the internal 
shock contributes the variable power-law spectrum, and the external shock 
makes the long-lived power-law spectrum. The central core of a massive 
star or the merged compact star gravitationally collapses into a black hole 
or neutron star, which somehow launches a collimated jet. Dissipation via 
shocks internally or with nearby matter under the photosphere determines the 
baryon loads and subsequent fireball evolution, where low baryon loads lead 
to a VHLF, as in Fig.[3]and TableQ 

factor Tc at = 
change in the baryon loads leads to a large difference in 
the Lorentz factor ATc/Tc > 10. 



lO'^r.s^''^ in Fig. [3] where a slight 



We note that the bulk Compton emission is another possibil- 
ity to produce the extra GeV component in the VHLF models. 
Since the emission depends on the ambient photon fields, we 
leave it as a future work. 

4.2. GeV spectrum: internal shock synchrotron model 

A VHLF could be a key for internal shocks to make the 
variable power-law spectrum that is slightly rising as vFi, oc 
i/O i_jyOA^ sometimes over ^ 7 energy digits from > 10 GeV 
down to < 20 keV. In the following, we propose the simple 
internal shock synchrotron model with a VHLF for the extra 
spectral component. 

The internal shocks convert the kinetic energy into internal 
energy with the energy density of 



U' = 



47rricr2 167rc3f,2r6' 



(54) 



where the shock at a radius r^h = IcT^ty makes a variability 
of time f,.. We assume that electrons are accelerated in the 
internal shock to a power-law distribution of Lorentz factor 
7g, dn'^/d^[, oc 7^ for 7^ > 7/,, and p > 2. As we have dis- 
cussed in § ^3.3| and |3.5l abundant pairs likely exist. We 
parametrize the e^ number density as 



: Tin' = n- 



U' 



ympC^ 



(55) 



where 7' is the random Lorentz factor of protons. If a fraction 
of internal energy goes into the electron acceleration, = 
eeU', the minimum electron Lorentz factor is given by 



p-2mp 7' 
p-1 m^TZ 



300 e,j'n-'fp, 



(56) 



where fp = 6{p-2)/{p-l). We further assume that a fraction 
of internal energy goes into the magnetic field amplification. 



B' = (SttcbU' 



1/2 



3,1/21/2 ,-1 



(57) 



Thus, the electron synchrotron cooling is effective above 



UTB'^-Vt,, 



-4x10^ VlnW-^et 



"=5,-2 ■ 



(58) 



With the synchrotron formula (^(7^) = 3qeB'Tj'^ /AirmeC, we 
have the synchrotron characteristic and cooling frequencies 



v„=5Q eV TfL%\}_,elef_^i^n-^fl, 



3/2, -3/2 
'I'.-Stfi _2 I 



(59) 
(60) 



respectively. Because of the slow cooling < Vc for a VHLF, 
the spectrum is 




vFT oc 



where the synchrotron luminosity a.\.v = Vc is about 



^syn(^c) 



p-1 



(p-2)/2 



(61) 



(62) 



since electrons above 7^ cool effectively. 

Figure |5] shows the internal shock synchrotron spectrum 
with a VHLF, F = 10'*, for fiducial parameters, L= 10^^ erg 
s-i, e, = 1, Cb = 10-2, f, = 10-3 s, 7^ = 10, 7' = 10, p = 2.2, and 
redshift z = 1 . Here, the attenuation by the creation with 
the cosmic infrared background becomes important at > 20 
GeV, which is taken into account using the best-fit model of 
Kneiske et al. (2004) (Kneiske et al. 2004). 

From Fig.|5] we can find that the VHLF internal shock syn- 
chrotron model has the following advantages. 

(1) First, the extra component has a single rising power- 
law spectrum over many energy digits. Thanks to the 
strong dependence of Vc on the Lorentz factor, oc F^, 
a VHLF internal shock has a high cooling frequency Vc 
beyond 10 GeV. The F dependence is oc F-^ for 
the characteristic frequency v,„, extending a power-law 
below < 1 keV. In combination, the rising segment of 
the vFjj oc jy(3-/')/2 spectrum is stretched to 



'2x 10' 



F 

104 



10 



(63) 



for fiducial parameters with Eqs. ( |59l ) and ( |60] l. There- 
fore, a VHLF could be crucial for making a single 
power-law spectrum that is rising over >7 energy dig- 
its. 

(2) Second, the VHLF models may have the high internal 
shock efficiency, which can supply sufficient energy to 
the extra component that is comparable to the photo- 
spheric Band component. From Eq. ( |62] |. the lumi- 
nosity of the high-energy component is almost equal to 
the photospheric one, Lsyn(i'£-) ^ i^ph ^ L, if p 2 and 
w 1 . The high-electron-energy fraction ~ 1 may 
be realized by highly efficient internal shocks between 
VHLF and slow shells, and also by e^-ri ch s hell s ex - 
pected for fireballs after dissipation (see § 33.3l and l3.5l ). 

Alternatively, the power-law spectrum may be shaped by 
the superposition of emission from continuous regions and/or 
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Fig . 5 . — Spectrum for the internal shock synchrotron model with a VHLF 
r = lO'*, for fiducial parameters, L = lO" erg s"' , = 1, = 10"^^ = jq-s 
s, 7?. = 10, 7' = 10, p = 2.2, and redshift z = 1, with (thick solid line) and 
without (thin s olid line) the attenuation by the cosmic infrared background 
IKneiske et alI 2004). Thanks to the dependences on the Lorentz factor of the 
cooling frequency Vc oc F** and of the characteristic synchrotron frequency 
h'ni oc F"^, the rising segment of the h'Fv oc jy<^"''V2 spectrum is stretched 
over > 7 energy digits. The maximum synchrotron cutoff is determined by 
the cool ing time t*™"' oc F in Eq. ill) or by the dynamical time uf^".^ oc F"'' in 
Eq. j73), which could produce the ~GeV cutoff as observed in GRB 090926 
for a VHLF case, F ~ 10*-10^. The spectral integrati on o ver continuous re- 
gions and/or times would lead to a flat spectrum in Eq. i64\ (dashed line). We 
also show the photospheric spectrum (solid line) with a Comptonized com- 
ponent (dotted line), and the upscattered photospheric spectrum with (thick 
dot-dashed line) and without (thin dot-dashed line) the attenuation. These 
are compared with the CTA sensitivity where we raise the public one by (50 
hr/10 s) for simplicity (as the worst case). The actual sensitivity in low energy 
could be worse than the public one l iTeshimall20 1 01) . 

times. Since the dependence of the cooling frequency i/c on 
the Lorentz factor F is strong, a slight change in F results in a 
large shift of I'c, and hence, in a relatively flat vF^ spectrum. 
As an example, we assume LccT^ and 7 cx F'/^ for simplic- 
ity. Then, with u, cx F^^L'^/^ and L,yn(iy,) cx L(L^T-^)p-^ in 
Eqs. (|59] l, (|60] |. and (l62T i. the envelope of the spectrum in- 
tegrated over continuous regions and/or times becomes rela- 
tively flat as 



2A+(A-9)(p-2) 
CX V 1<>-3A 



,0.05 



for ly < I'c, (64) 



where the last equality is for p = 2.2, and A = 1, i.e., the 
case when the luminosity is proportional to the kinetic energy, 
which may be reasonable. 

Our model predicts a temporal correlation between the ex- 
tra power-law component in the low (~keV)-energy region 
and the high (^GeV)-energy region, because a single emis- 
sion mechanism produces the whole power-law component. 
In the integrated synchrotron case in Eq. ( l64b . the low-energy 
component could be delayed by the variability time, i.e., the 
dynamical time for the Lorentz factor to change. The pho- 
tospheric Band component is also expected to be temporally 
correlated with the power-law component, at least within the 
variability time, because the photospheric luminosity is com- 
parable to the kinetic one for the VHLF range, as shown in 
Table [1] and the internal shock emission is delayed by the 
variability time at most; 

Toma et al. (2010) dToma et al.ll2010l) showed that the pho- 
tospheric emission can be efficiently scattered by the elec- 
trons in the internal shocks, and the Compton upscattered 
photospheric emission can explain the observed distinct high- 



energy component. The low-energy part of the power-law 
component at < 20 keV is attributed to the synchrotron emis- 
sion, i.e., a different origin from the high-energy part. This 
is different from the VHLF models that employ only a single 
emission mechanism (synchrotron emission) and may be fa- 
vorable for explaining a single power-law component over >7 
energy digits. In the VHLF models, the upscattered emission 
is beyond the Fermi energy range, because the characteristic 

frequency of upscattered photons is ^ 7'm£peak > 1 TeV (see 
Fig. |5]), which is rather a nice target for the future Cherenkov 
Telescope Array (CTA) Project. The electron cooling by the 
Compton upscattering is also suppressed by the Klein-Nishina 
effect. The emissivity ratio of the upscattered photosphere to 
the synchrotron is 



Pupil'e) Kb < ''KNh'j] 



(65) 



where fKN(7c) = rsh'«cC^/7e and the step function approxima- 
tion in the second equality is appropriate if d\vJJ '^^j d\s\v < 2 
for J/ > i^KN dWang e t al. 2010; Nakar et al. 2009). Since the 
comoving energy density of photospheric photons is approxi- 
mated by 



f/^h < ^^kn(7:)] ^ U' 



^'KN(7e) 



-peak 



2-Q 



(66) 



for i^KN < £peak. the upscattering dominates the synchrotron in 
the electron cooling yup(7F) > 1 at 7^ < 7^ where 



7f' 



TnigC^ -l/(2-a) J. in5 -n -l/(2-a) .^-s 

-e^ --5x 10^r4epeak.MeveB:2 (67) 



£peak 



for a = 1 . This is lower than 7c in Eq. (158b . and hence, the 
electrons mainly cool via synchrotron, not via the inverse 
Compton, for the VHLF models. (Note that we have to take 
the Compton cooling into account when we precisely calcu- 
late the low-energy end of the integrated spectrum in Eq. ( |64l ), 
which is produced by the d ecelerated shocks.) We can also 
negle c t flie SSC emissi on dPilla &Loebl 119981: ISari &EsinI 
120011; iFanetallEooi; iZou et alJ l2009riNakar et all 120091: 
Wang et al. 20iy), since the synchrotron photon density is 
lower and the peak energy is higher than those of the pho- 
tospheric emission. 

The synchrotron component could also be detected in the 
optical band. It could explain the bright optical prompt emis- 
sion in some GRBs, w hich is brighter than the extrapolation of 
the Band component dYost et alfcoOTi; iRacusin et alj|2008l) . 

4.3. GRB 090926: creation cutoff, cooling break or 
maximum synchrotron cutojf? 

Recently, a spectral break at approximately 1.4 GeV has 
been found in the extra power-law component of GRB 
090926. We consider the possible origins of the break in the 
following. 

(1) The 77 annihilation with e^ creation leads to a spec- 
tral break ( Baring & Harding 1997; Lithwick & Sari 200ft 
iRazzaque et al.ii2004t iMurase & lokaii2008i: lAoi et abi2009[) 
at 



u^-f 20 GeV 



._1 1 2+2/3 



2-/3 

peak.MeV 



(/3-1) 



,(68) 



which is determined by the optical depth condition 
t{v) ~ imn'Jv' > v')aTr,bl2T = 1, where v' = mYjv', 
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- 1(0 -1)/\6B ^^\0 + 1)] - 0.1, (iSvenssonl fT987t 
iMurase & Iokall2008h and we approximate the photon num- 
ber density as 



«>;>z>') = 



(69) 



with /? = 2.5. Thus, the break at ~ 1 .4 GeV suggests the bulk 
Lorentz factor of ^ 600. This constraint is applicable to the 
outermost shell, not excluding VHLF shells behind the slower 
shell. (Note that this alignment is eve n pre ferred by the con- 
siderations of the GeV onset delay in ^4.4l i. The creation 
cutoff disappears if D' < £peak, that is. 



167rc2f,,£peak(/3-l) 



1/4 



'3 X 10' L't 



3 ,1/4,-1/4-1/4 



53 'i'.-3 ^peak.MeV 



(70) 



in the VHLF range, because the photon density is almost con- 
stant below Epeak for the usual spectral index. We note that the 
exponential cutoff is usually sm oothed to a broken po wer-law 
by multizo ne effects jAoi et al. 2009; Li 2010a; Gran ot et alj 
|2008; Bos niak et alj|2009l) . 

(2) The second possibility is the cooling break at 



ly, = 9 GeV Tll/^^^t, 



,g-3/2 
-3^^,-2 ' 



(71) 



in Eq. (|60l). Thus, the break at 1.4 GeV suggests a VHLF 
of r ^ 10"*. The change in the spectral index at the cooling 
break is 0.5 in Eq. dHTT i. which can be used as a test of this 
possibility, although the photon number is insufficient in the 
current observations. If the low-energy part below the cooling 
break is produced by the spectral integration in Eq. ( |64l ), the 
change of the index is smaller than 0.5. 

(3) The third possibility is the maximum synchrotron cut- 
off. 



3q,B' 



cool / -'^e" I 
V,^.,^ = CPs- 7„ 



27 m^c r 
167r K 



= 50 K-T3 GeF,2) 



which is only dependent on the bulk Lorentz factor F, and is 
determined by the balance between the acceleration time t!^^^ = 
K^'^nigC / q^B' and the cooling time t^^^^ = 3m(.c/4(T7-f/^7^', so 

1 /2 

7max = (6TTqe/ KajB'^ , whcrc the coeffi cient (ps = 0.2294 
is quo ted from Wijers and Galama (1999) jWiiers & Galamal 
[T999h . If the break at - 1 .4 GeV is due to the maximum syn- 
chrotron cutoff limited by the cooling time t',™"', the shock 
acceleration has to be much slower than the Bohm limit, 
K ~ 50 r3 ^ 1, i.e., the scattering mean free path is much 
larger than the Larmor radius. Apart from the break in GRB 
090926, the VHLF models (F - 10'*) predict the maximum 
synchrotron cutoff in the TeV region for k ^ 1, which is a 
nice target for the future Cherenkov Telescope Array (CTA) 
Project (see Fig.|5]l. 

On the other hand, if < Vc, i.e., F > 10"*, the max- 
imum synchrotron cutoff is limited by the dynamical time 
Ff,, > f^^^, rather than by the cooling time, yielding 



'dyn 



7max = qeB'Tty/ nnieC, and hence. 



,,dyn _(-) 1 (-< \j p-6r3/2 3/2 -3 -2 



(73) 



Therefore, the break at ~ 1.4 GeV could be produced by 
VHLF flows with F ~ 6 x lO'*. It is a unique feature for the 



VHLF models to be able to accompany the maximum syn- 
chrotron cutoff limited by the dynamical time, not by the cool- 
ing time. 

(4) The last possibility is that the extra component of GRB 
090926 might be the annihilation line from the photo- 
sphere at the blue shifted energy. 



v± = TnieC^ - 0.5 GeV F3 



(74) 



which is broadened by the order-of-unit y distribut i on of 
the Lorentz fa ctor on the photosphere jloka et al.l 120071; 
IMurase &Iokall2008 ). This scenario might be possible if 
the pairs are continuously created on the photosphere, al- 
though the mechanism for the creation is not apparent. 

4.4. GeV onset delay 

Fermi discovered that the high-energy emission (> 100 
MeV) is delayed behind the onset of the MeV emission in 
almost all LAT GRBs. The delay time in the rest frame 
is fdeiay ^ 1 s for long GRBs and 0.1 s for short bursts, 
GRB 081024B and GRB 090510. These delays are not just 
caused by the flux increases above the LAT detection thresh- 
old, but by the spectral changes in the Band and/or extra com- 
ponents at least in the well-observed bursts. Since the ob- 
served delays of the extra component are larger than the vari- 
ability timescale of the Band component, ^ 0.01-0.1 s, the 
physical origin of the delay is not likely the kinematic effect 
(lTomaetal.il2010h . 

In the hot photosphere-internal-external shock model in 
^4.1l and Fig.|4] the GeV delayed phase arises when the emis- 
sion from internal and external shocks is weaker than the pho- 
tospheric emission, i.e., almost all the energy escapes in the 
form of the photospheric luminosity, not the kinetic luminos- 
ity. This is realized for the dimensionless entropy in the mod- 
erate range ry, ~ 10^ < V < Vki ^ 10'* according to the dissi- 
pative hot photosphere model in Table[T] Fig. [3] and Eqs. ( l23T l 
and ( |36] |. On the other hand, in the GeV bright phase, the ki- 
netic luminosity is comparable to the photospheric luminos- 
ity, i.e., in the high-entropy range ry^i ^ 10"* < 77 < ?7<:3 ~ 10^ 
in Table [T] Fig. [3] and Eqs. (|36] | and (l38T l, which is also the 
VHLF range. 

Therefore, in our picture, the baryon loads decrease pro- 
gressively from the GeV delayed phase to the GeV bright 
phase, across the critical entropy rj ^ ijki in Eq- (l36] l and 
Fig-E] As argued in § ^3.1l and l3.2l the baryon is entrained at 
the dissipation radius r„„ which is a function of 77 in Eq. (|35] |. 
Then, we can predict the delay time as the light crossing time 
of the dissipation radius for 77 = ?7/n. 



^delay ' 



rmiVki) rt 



'05? A^/^T"^/^ r*/^ ('75') 

U.J A ^eOOkeV^ s 



where is the base size of the flow in Eq. (|9]). Interestingly, 
this predicted timescale is comparable to the observed de- 
lay time fdeiay ^ 0.1-1 s. This coincidence also supports our 
picture that the VHLF fireballs are responsible for the extra 
high-energy emission. In addition, Eq. ( fTST i combined with 
the Epeak-L Yonetoku relation in Eq. ([T]) gives 



^delay ' 



niVkl) 



oim) 



0.5sL.\'^ry^, (76) 



which has a weak dependence on the luminosity, consistent 
with the observations. Further observations of the delay time 
would discriminate models since other models have differ- 
ent parameter dependences. For example, if the Coulomb 
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collisions control the kinetic luminosity, the delay time is 
fdeiay ^ r,„{rjc) oc L'^'l^v]^'^ f]l^ with the critical entropy 77c in 
Eq. ( fSTT ). The delay caused by the creation cutoff also has 
a different dependence fdeiay « LT~^ (lLill2010ah . If we can re- 
fine the delay time measurements in the future, the delay time 
might be used as a distance indicator like the Epeak-i' Yonetoku 
relation. 

The delay timescale is also comparable to the light crossing 
time of the progenitor star 



fdel; 



ay 



star 
C 



'03 s 



10'" cm 



(77) 



Therefore, a natural picture is that the fireball dissipation via 
baryon loads is controlled by the environment just outside the 
star As the environment changes with time ^ /^stai/c, the 
baryon loads decrease and the dissipation radius increases be- 
yond the thermalization radius, leading to a VHLF fireball via 
the coUisionless bulk acceleration in ^3.41 and hence, to the 
high-energy emission with the onset delay. 

Although we have used the light crossing time to estimate 
the timescales, the actual environment could be more com- 
plex, which may not be governed by the light speed. After the 
jet breakout, the jet is likely surrounded by a cocoon, which 
consists of the decelerated jet and the shocked stellar enve- 
lope. The velocity of the cocoon depends on the jet luminosity 
and the stellar structure, although it is typically Vc/c^O.1-1 
(fToma et al. 2007). If the accelerated particles are important 
for the baryon loading, the typical velocity is the light speed. 
The baryon loading process is n ot cl ear, and further discus- 
sions are given in § ^5.3ll54l and |5.5l 

5. PREDICTIONS AND OPEN ISSUES 

5.1. TeV neutrino 

We can predict ^TeV neutrinos and their temporal an- 
ticorrelation with extra GeV high-energy 7 -rays in the 
hot photosphere-internal-external shock model in 34.11 and 
Fig. |4] In this model, the baryon (protons) is loaded 
below the photosphere, leading to the fireball dissipa- 
tion as suggested by observations (see 33.11 1. The en- 
trained protons are relativistic in the comoving frame of 
the shocked fireball with the Lorentz factor 7' ^ r,„/2r.5 ^ 



^Jr]/Vsl2 in Eqs. (l20l l and (l35l l. The pp collisions be- 
tween these pro tons produce pions, wh i ch immediately decay 
into neutrinos d Paczvriski & Xu 1994; Derishev et al) 119991: 
iBahcall & Miszaros 2000; Meszaros & Rees 2000b) via 



7r~ -> /i~ + e~ + T>e + + . 



(78) 
(79) 



Each neutrino shares ^ niT^/Amp ^ 5% of the primary proton 
energy, so that the observed neutrino energy is 



4m „ 



. mlpnipC 



'0.2 TeV (^) 



(80) 



The neutrino luminosity originates from the proton kinetic lu- 
minosity, which also produces the extra GeV 7-ray compo- 
nent in our model. That is, the same kinetic energy is shared 
by neutrinos and extra 7-rays. Therefore, the ~TeV neutri- 
nos are predicted to anticorrelate with GeV 7-rays, as shown 
in Fig. |6] and Table [T] Such a temporal prediction would be 
interesting for the upcoming multimessenger astronomy. The 
neutrino fluence may be comparable to the Me V 7-ray fluence 



since the kinetic energy is comparable to the photospheric en- 
ergy, as shown in Eq. ( fTsT i and Table [T] However, as dis- 
cussed in 34.41 the pp collision is effective only in the GeV 
delay phase with a timescale ~ fjeiay 1 s less than the to- 
tal duration T ^ 20 s. Therefore, the neutrino fluence would 
be ^ fdeiay/r ^ 0.05 timcs less than the MeV 7-ray fluence 

J"- 10"^ergcm-2. With the GRB event rate T^grb 10^ yr"', 
the diffuse neutrino background from GRBs is estimated as 



1 ^delay ^ ^ 

A-K T 



X 10 



r'°GeVcm-^ s"' sr''^!) 



which is less than th e current limits dAbbasi et alj20ldll2009t 
iFukuda et al.lll992h and less than the IceCube design sensitiv- 
ity by an order of unity. If the GeV bright bursts are a minor 
population as the LAT bursts are ^ 7% of Fermi bursts, the 
diffuse neutrino flux is higher by ^ T/fdeiay ^ 20, i.e., com- 
parable to the diffuse 7-ray flux from GRBs, which could be 
detectable using IceCube in the near future. We note that Ice- 
Cube is sensitive to neutrinos above TeV energy, and not so 
much to sub-TeV neutrinos in Eq. dSOl l. However, a fraction 
of protons would have the random Lorentz factor larger than 
7^, in Eq. dSOl l. as discussed in 33.21 so that the detection is not 
completely hopeless. 

The photomeson interactions (777 — ^ mi) could also gener- 
ate pions, and hence, neutrinos if the photon spectrum under 
the photosphere is already nonthermal, as discussed in 33.51 
Each neutrino energy is also similar to the pp collisional case 
in Eq. dSOb , since the primary proton energies are the same. 
We also expect the anticorrelation between ~'TeV neutrinos 
and GeV 7-rays in Fig. |6] as in the pp collisional case be- 
cause the energy source is again the proton kinetic energy, 
which also produces the extra GeV 7-ray component. How- 
ever, the neutrino luminosity is suppressed by the factor in 
Eq. ( |49] l compared with the proton kinetic luminosity because 
the Bethe-Heitler process also occurs simultaneously and con- 
sumes the proton energy into creation. Therefore, the pho- 
tomeson neutrinos may be difficult to detect unless the photon 
index is hard /? < 2 in Eq. ( l49b . 

The neutrino emission could precede the 7-ray emission 
when the jet is sti ll inside the progenitor star jRazzaque et alj 
I200ll^ueil200l . In this case, the accelerated protons inter- 
act with matter of the progenitor star or synchrotron photons. 
The 7-rays cannot escape owing to the large optical depth. 
The time delay between neutrinos and 7-rays is expected to 
be about Rstm-/vj ^ 1-10 s, approximatel y the time taken b y 
the jet to emerge f rom the progenitor s tar jToma et al. l2007h. 

Wa ng and Dai jWang & Pail 120091) and Murase jMurasd 
I2OO8I) also discussed the high-energy neutrino emission from 
the dissipative photospheres of GRBs. We can predict that the 
high-energy neutrinos in their models would also temporally 
anticorrelate with the GeV 7-rays. 

5.2. MeV "f-my spectrum 

In the photosphere model, it is unclear how to produce the 
high- energy nonthermal tail of the Band spectrum (loka et al] 
l2007h . The main problem is the source of the nonthermal 
energy that is comparable to the total energy. The nonthermal 
tail has to be produced near the photosphere ~ 1, since the 
spectrum is thermalized above the observed frequency. 
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Fig. 6. — The anticoiTelation between ~TeV neutrinos and extra variable 
GeV 7-rays is schematically shown. This is predicted independently of the 
neutrino generation processes, either pp or photomeson interactions, because 
the same proton kinetic energy is shared by neutrinos and extr a GeV 7-rays 
in the hot phot osph ere-internal-extemal shock model in j|4. H and Fig.|4] In 
our picture of j|4.4l the GeV onset delay phase corresponds to the baryon-rich 
(low entropy) phase, associated with the effective p p ne utrino emission with 
~TeV energy in Eq. (80) and the diffuse flux in Eq. 



through Compton scatterings by the thermal electrons with 
temperature less than that of the nonthermal photons. This 
is because a photon changes its energy by Ai''^ jv'^ ^ - 

v'^)/meC^ in a single scattering, so that a relatively large opti- 
cal depth Tj ^ v'^j Av'^ ^ nieC^ /v'^ is necessary for thermal- 
ization at v'^ > AkT^. Note that the necessary optical depth 
is not tt ^ (v'^l /S.v'^)^l^ here, since the fireball is expanding 
with decreasing tj. 

As a by-product of the discussions in the previous sections, 
we could find a hint for the energy source of the high-energy 
nonthermal Band spectrum in the photosphere model. That is 
the relativistic baryon (protons) component loaded below the 
photosphere at th e fireball dissipation, as suggested by obser- 
vations (see ^3. 11 1. The entrained protons are relativistic in the 
comoving frame of the shocked fireball (see ^3.21 ). and hence, 
can heat electrons via pp collisions (see 33.31 ). p7 (Bethe- 
Heitler and photomeson) processes, and Coulo mb c ollisions 
with created by pp and p"f processes (see 33.5b . Then, 
the heated can produce the nonthermal broken power-law 
spectrum through the unsaturated Comptonization of the ther- 
mal photons. The total proton energy before thermalization 
is just comparable to the photospheric radiation energy in 
Eq. ( fTSl l, as required by the nonthermal tail of the Band spec- 
trum. This energy equipartition comes from the energy and 
momentum conservation in Eqs. (fTTl i and (fTSl l without fine 
tuning. For the electron heating near the photosphere tj ^ \, 
the VHLF range 77 > 77^1 may be preferred because the dissi- 
pation radius r,„ is near or above the baryonic photosphere rph 
(see Fig.©. 

Alternatively, the plasma turbulence could be initiated at 
the baryon loading, at least with mildly relativistic velocities. 
Since each fluid element has a relative velocity, the turbulence 
would scatter thermal photons to a broken power-law spec- 
trum, where the F-parameter can be of order unity for the dis- 
sipation near the photosphere. The turbulence will be damped 
within the eddy turnover time. However, once a nonthermal 



spectrum is formed, the spectral shape is almost preserved be- 
low I'th in Eq. (l82T i. 

Beloboro dov (2 09) (iBeloborodovl 120091) (see also 
Ref. Derish ev et al.l (119991) ') has recently shown that the pn, 
pp, and Coulomb collisions can heat to produce the 
nonthermal Band spectrum in the context of the neutron- 
loaded fireballs. Since the neutron models only produce 
mildly relativistic nucleons, it is interesting to calculate the 
case of the relativistic protons and neutrons in the shocked 
fireball frame. In contrast to the neutron models, the p^ 
(Bethe-Heitler and photomeson) processes are also important 
for the relativistic case as shown in 33.51 The 777 processes 
have an advantage that the random Lorentz facto r of protons 
does not drop to unity in most cases (see 33.51 ), and hence, 
the fireball can realize both the nonthermal photospheric 
spectrum and the fair fraction of kinetic energy for the 
internal and external shocks at the sa me time, in contrast 
to the pp collisions. Kazana s et al. (iKazanas et a l.l 120021: 
iMastichiadis & Kazanasl20^9b also considered the relativistic 
proton accumulation and the Bethe-Heitler process in the 
supercritical pile model, which may be relevant to the late 
evolution above the photosphere. 

5.3. Origin of the £peak-i Yonetoku relation 

Thus far, we have just used the Epeak-^ Yonetoku relation 
in Eq. ([T]) as an empirical relation. Let us consider the phys- 
ical origin of this relatio n in the dissipative hot photosphere 
model. As discussed in 33.11 the Cpeak-i' Yonetoku relation is 
reproduced by the fireball dissipation under the photosphere 
that is associated with the deceleration of the fireball, proba- 
bly via the baryon loading. The necessary (isotropic) baryon 
loading rate M can be obtained from Eqs. (|9]), fT% . and ( 1351 ), 
without using the Yonetoku relation, as 



M' 



LT, r 



10"' Mq s 



r2 j-A 

m. 10^53 ^eOOkeV- 



(83) 



Interestingly, the dependence L^T~* disappears if the epeak-^ 
Yonetoku relation in Eq. ([T]) is satisfied. The remaining de- 
pendence in Eq. ( [83] ) is the Lorentz factor of the slow mass F, 
and the baryon loading radius (dissipation radius) r„,, which 
are all determined by the environment that supplies the baryon 
into the jet. In other words, the Epeak-i' Yonetoku relation is 
satisfied if the baryon loads are controlled by the environment 
and do not differ event by event so much. This picture is also 
consistent with the fact that the GeV onset delay is controlled 
by the environment in 34.41 However, it is still difficult to 
derive the normalization of the baryon loading rate, partly be- 
cause the current MHD numerical simulations have not yet 
implemented the key physical processes of radiation trans- 
fer and creation dLazzati et alj 120091: IZhang et all 120041: 
iMizuta et ai]|2004 120101) . For Eq. ([831), the essential relation 
is r,„ oc y/rj in Eqs. ([TTb and ([35] ). which states that the fire- 
balls are relativistic even after the dissipation, derived from 
the new Equations ( [TT] ) and SlTi . 

5.4. Baryon load and jet structure 

As discussed in 33.11 the epeak-i- Yonetoku relation in 
Eq. ((Tj) suggests the fireball dissipation under the photosphere 
that accompanies the deceleration of the fireball, probably 
via the baryon loading. However, the actual baryon loading 
process is unclear, partly because the current MHD numer- 
ical simulations have not yet implemented the key physical 
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processes of radiation transfer and creation ([Lazzati et alj 
I2OO9I: iZhang etaPllOOl IMizuta et al.ll2004l IIoIoIk 

From the kinematical viewpoints, we may argue that the 
baryon loading process would require the nonspherical con- 
figuration. Otherwise, the radiation-dominated outflow from 
the central engine with the initial width ro ~' 10^ cm cannot 
be fully dissipated at the dissipation radius r,„ in most cases. 
We can show this by comparing r,„ with the radius where the 
rapid outflow with width ro ^10^ cm completely runs into the 
shocked region with the Lorentz factor r„, as 



rb/ro 



(84) 



with Eqs. ^ and ( |35] ). This ratio is larger than unity, i.e., the 
dissipation is not completed at r,„, for almost all parameter re- 
gions -q > lO^r^' in Fig. |3] This problem of the incomplete 
dissipation applies to all dissipative photosphere models with 
r > 10^ and Lph - L. In other words, the one-dimensional 
picture of the spherical reverse and forward shocks is insuf- 
ficient to explain the GRB observations in the photosphere 
model. The picture should be rather close to the complete 
merger with turbulence. 

The nonspherical configuration may be plausible because 
the low-entropy fireball is pushed by the high-entropy fire- 
ball and the system is subject to th e Rayleigh-Taylor insta- 
bility and/ or the plasma instabilities (" Waxman & Piranlll994t 
lloka et ai] 12007). Then, the full dissipation condition (r^ > 
r^/o in the above) would be replaced by r„, > F^ro that 
is easily satisfied. One may also consider a jet configura- 
tion of the fireball that e ntrains baryon (protons and neu- 
trons jLevinson & Eichlei1i200 3)) through the boundary of the 
jet. If the central engine works intermittently, the baryon sur- 
rounding the jet, previously in pressure equilibrium with the 
jet, will enter the funnel sideways. As long as the jet opening 
angle 9j{> F^') is smaller than the causal angle F~', the fire- 
ball can be fully dissipated via turbulence, e. g., caused by the 
Kelvin-Hebnholtz instabil ity jO oodman & MacFadvenll2008l; 
ISironi & Goodmanll2007t IZhang et al.ii2009i) . Note that the 
jet opening angle remains constant after the jet reacceleration 
because the causal angle shrinks all the time. 

Although the causal angle F~' ^ 0.01 7^4 '''^Fj '''^ for the high 
entropy flow with r/ > 10'* may be smaller than the conven- 
tional opening angle dj 0.1, the jet could have two com- 
ponents or a continuous structure with the high entropy flow 
surrounded by the low to moderate entropy flow of 77 ^ 10^- 
10'*. We note that the two-component jet is suggested for 
GRB 080319B (Racusin et al. 2008). If this is the case, the 
extra GeV component associated with the high entropy flow 
(see Table ^ may be observed only from the viewing angle 
near the jet center, possibly consistent with the fact that not 
all observed bursts are bright with GeV 7-rays. 

The initial variability time ^ ro/c ^ 3 x 10""* s reflecting 
the central engine size could be preserved at the dissipation 
since the successive shells may not contact with each other 
for r,„ < roF^ in Eq. (|84] i. If the dissipation is nonspherical as 
discussed above, the timescale could instead be determined 
by the crossing time of the causal region, rm/cTm = rh/c ^ 
3 X lO""* s, with Eqs. dUl and 

5.5. Multiple baryon loads 

The baryon-loaded shell could expand and collide before 
coasting against another slow mass Mj2 with a Lorentz factor 



Fj2. This second merger is characterized by the dimensionless 
entropy 772 = £','F,./Mj2C^, since the rapid shell is radiation- 
dominated as in the first merger. The merger is divided into 
three types: 

( 1 ) If 772 > /Fj2, we may neglect the second merger since 
the rapid shell is not sufficiently decelerated and the 
fireball temperature is almost constant. 

(2) If ri2 < the total amount of baryon is determined by 
the second merger. The Lorentz factor of the fireball 
goes down to F„j2 ~ V^siVi and then up to Fc2 ^ 772 
with a different base radius rb2 ^ rni/^mi- We may 
consider that the initial condition of the dissipated fire- 
ball is reset by the second merger. 

(3) If 77 < 772 < F^/F52, the total amount of baryon is not 
changed so much, while the fireball is decelerated to 
have a different temperature and base radius. The ini- 
tial radius of the dissipated fireball is not r„j/rh ^ F„, ^ 

/TV/, but r,„2/rfc 



r,„2 ^ -yF.52772 that tends to bring 
the dotted line upward in Fig. [3] As a result, we have 
different critical entropies from 77^1 , 7/^2 , ila in Eqs. (l36t , 
(|37] |. (|38] ). and therefore, a different maximum Lorentz 



factor from Fc.max, r^:^,^^ in Eqs. (|4T]l, (gUi. It is 
straightforward to obtain these quantities in this case. 

We note that the last value of the base radius determines the 
observed relation in Eq. 

If the merger type is (3), the necessary baryon loading rate 
becomes different from that in Eq. (|83] |. Therefore, the merger 
type (3) could introduce an outlier in the Yonetoku relation, if 
the origin of the Yonetoku relation is correct in ^5.31 To have a 
tight Yonetoku relation, the baryon loading needs to be larger 
at the outermost radius, which may be reasonable since the jet 
boundary also becomes larger 

5.6. Early X-ray afterglow 

Swift discovered the s teep and shallow decay phase in the 
early X -ray afterglows (Zhang et al. '2006'; loka et al. '20061; 
Panaites cu et al. 2006; Zhang et al. 2006; Huang et al. 200g; 
Sato et alJ l2007). The important point of these observations 
is that the early X-ray afterglows are too dim for the con- 
ventional internal-external shock model, in which the kinetic 
energy left for the afterglow emission is comparable to or 
usually larger than the prompt energy released at the inter- 
nal shocks, so-called the internal shock efficiency problem 
(lokaetal. 2006). In the hot photosphere-internal-extemal 
shock model in Fig. H] and Table [T] this is not a problem be- 
cause the kinetic energy fraction of the afterglow may be rel- 
atively small for the moderate entropy range 7/*(^ 10^) < 77 < 
77^1 (^ 10"^) in Table [U where the fireball becomes optically 
thin in the accelerating phase before converting all the radia- 
tion energy into kinetic energy. 

On the other hand, the LAT bursts with bright GeV 7-rays 
do not seem to be associated with the steep and shallow X- 
ray afterglows, although more events are necessary to con- 
firm this fact. In these bursts, the kinetic energy fraction of 
the afterglow may be larger than that of the Swift bursts. This 
may be consistent with our picture that the extra GeV emis- 
sion is produced by the relatively high entropy fireball with 
V > Vkii^ 10'*) in the VHLF range, whose kinetic luminos- 
ity is comparable to the photospheric luminosity, as shown in 
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Table[T| If the jet is structured with the high entropy flow sur- 
rounded by the moderate entropy flow as discussed in 35.41 
the shallow decay in Swift events may be detected for ob- 
servers in the direction of the moderate entropy flow, since 
the high-entropy region with large kinetic energy is progres- 
sively seen as the afterglow is decelerated Qbma et al. 2006|). 

5.7. Short GRB: magnetized jetfmm white dwarf? 

Short LAT GRBs 081024B and 090510 have smaller GeV 
delay time, fdeiay ^0.1 s, than long GRBs. According to 34.41 
this implies that the size of the progenitor for short GRBs is 
smaller than that for long GRBs. Interestingly, the inferred 
size cfdeiay ^ 10^ cm is comparable to the white dwarf radius. 
The white dwarf model may also be favorable to explain the 
extended emission observed in the short GRBs, since the ac- 
cretion time of the white dwarf material into the central en- 
gine can be much longer than the neutron star case, where the 
neutron star accretion is too fast to ex plain the long timescale 
(> 100 s) of the exten ded emission (iBarfliebny etal J 120051: 
iNorris & Bonnelill2006h . 

For short GRBs, the LAT (~GeV) fluence is comparable to 
and even larger than the GBM (^MeV) fluence, in contrast 
to long GRBs with th e LAT-to-GBM fluence ratio of 0.1 
jGhisellini et all 120101) . This may require a modification to 
the simple version of our model, because the extra GeV lu- 
minosity is usually less than the photospheric luminosity in 
Eq. (|62] |. One possibility is that the short GRBs could have a 
low entropy phase with rj < 77, (baryon-rich phase), in which 
the photospheric emission is suppressed below the kinetic one 
(see Table[T]i. The other possibility is a magnetized fireball jet 
because magnetic fields are not radiated away at the photo- 
sphere, reducing the photospheric emission (see 33.7l i. The 
magnetic fields also increase the kinetic energy by pushing 
the matter even above the photosphere via the magnetic pres- 
sure, which could enhance the GeV emission. 

6. SUMMARY 

We have investigated the fireball model with dissipation un- 
der the photosphere, as suggested by the observed spectral 
£pe-dk-L Yonetoku relation in Eq. ([T]i ( 33. 11 1. We find that the 
fireball can entrain the relativistic baryon component at the 
dissipation, with energy comparable to the radiation compo- 
nent, as derived from the energy and momentum conservation 
in Eqs. (fTTl i and (fT2] i. The relativistic baryon (proton) compo- 
nent can alter the GRB fireball dynamics and spectra in novel 
ways, which we have studied in this paper. Our main results 
may be summarized as follows. 

• As outlined in ^ the relativistic baryon component 
can reexpand to a very high Lorentz factor (VHLF) 
r ~ 10^-10^, much larger than the conventional upper 
limit r < 10^ (El Figs. [1] and O. Since the pres- 
sure is provided by the relativistic collisionless motions 
of protons (and the magnetic field generated by these 
protons), this mechanism may be called the collision- 
less bulk acceleration. The VHLF is achieved when the 
baryon load is low (i.e., high entropy) without efficient 
thermalization of protons via pp, p^ (Bethe-Heitler and 
photomeson). Coulomb, and plasma interactions (§ 33.31 
and 13.5b . The kinetic energy can become comparable 
to the total energy in the VHLF fireballs (EH Table[T] 
Fig- Ell, which can alter the internal and external shock 
emission (see below). These VHLF models are cur- 
rently consistent with previous observations ( 33.6b . 



• The VHLF fireballs can explain the extra variable GeV 
component in the GRB spectrum using the simple inter- 
nal shock synchrotron emission. In the VHLF models, a 
single emission mechanism (synchrotron emission) can 
make a rising power-law spectrum over >7 energy dig- 
its, and the internal shock efficiency can also be suffi- 
ciently high to produce the emission energy comparable 
to the total one, as observed by Fermi/LAT. The VHLF 
emission is also a nice target for the future Cherenkov 
Telescope Array (CTA), since the creation cutoff 
goes beyond the TeV range in the VHLF models. The 
main Band component may be attributed to the pho- 
tospheric emission, and the long-lived GeV compo- 
nent to the external shock in this hot photosphere- 
internal-external shock synchrotron model (§ 34. H and 

10 Fig, m. 

• In the VHLF models, the spectral break at ^ 1 .4 GeV 
observed in the extra component of GRB 090926 i s no t 
caused by the creation cutoff with F ^ 600 ( 34.3b . 
The spectral break could be the synchrotron cooling 
break for T ^ 10"^, or the maximum synchrotron cut- 
off, particularly limited by the dynamical time (not by 
the cooling time) for T ^ 10^. 

• The observed GeV onset delay can be naturally ex- 
plained in the dissipative hot photosphere model be- 
cause the baryon loads at the dissipation would be rich 
shortly after the jet breakout, and hence, the relativis- 
tic baryon component is thermalized into radiation in 
the fireball, leaving little kinetic energy for the GeV 
emission ( 34.4b . The timescale of the GeV onset de- 
lay can be calculated from the theory as the timescale 
at the dissipation radius, i.e., the timescale to change 
the baryon loading rate, and hence, the environment 
around the progenitor star The predicted delay time 
^ 0.5 s and also its weak dependence on the luminosity 
fdeiay « L"'/^ ^ L"'/^ in Eq. (|76] | are consistent with the 
observations. 

• The dissipative hot photosphere model predicts ^TeV 
neutrinos in Eq. (ISOl l and the anticorrelation between 
^TeV neutrinos and the extra variable GeV 7-rays, in- 
dependently of the neutrino generation processes, ei- 
ther pp or photomeson interactions, because the same 
proton kinetic energy is shared by neutrinos and extra 
GeV 7-rays. In the optimistic case, the diffuse neutrino 
background from GRBs in Eq. (ISTb could be detected 
using IceCube in the near future. 

• The relativistic baryon component loaded into the 
fireball at the dissipation could work as the as-yet- 
unknown energy source to deform the photospheric 
thermal spectrum into the nonthermal Band spectrum 
( 35.2b . Without fine tuning, the relativistic baryon com- 
ponent has the right amount of energy (compara ble to 
the thermal energy) to make the nonthermal tail ( 33.2b . 
The protons could heat via pp, pj (Bethe-Heitler 
and photomeson), and Coulomb interactions, leading to 
the Comptonization of the thermal photons into the bro- 
ken power-law spectrum. Alternatively, the mildly rel- 
ativistic plasma turbulence induced by the baryon load- 
ing could make the nonthermal spectrum. Further in- 
vestigations in this direction appear interesting. 
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• The spectral Epeak-i' Yonetoku relation can be repro- 
duced by the fireball dissipation if the baryon loading 
rate is determined by the environment that is nearly 
identical to any bursts in Eq. ( |83] ). Nearly identical en- 
vironments could be plausible if GRBs are produced 
by a certain class of progenitor For the derivation of 
the £peak-i Yonctoku relation, it is important that the 
fireballs are relativistic even after the dissipation, as de- 
duced from the energy and momentum conservation in 
Eqs. O and (fT2T i. 

• The actual baryon loading process at the dissipation has 
not been clarified. From the kinematics, the dissipation 
process would require the nonspherical configuration, 
such as the turbulence and/or the jet boundary ( ^5.4b . 
In the jet case, the causality argument suggests a two- 
component or structured jet. 

• A structured jet with both the radiation-dominated and 
VHLF flows could explain the steep/shallow decay of 



early X-ray afterglows observed in the Swift bursts (as 
viewed from the radiation-dominated flow) and the pos- 
sible paucity of the steep/shallow decay phase in the 
Fermi/LAT bursts (as viewed from the VHLF flow) in a 
unified fashion ( ^5.6t . 

• We have speculated that the short GRBs might accom- 
pany the magnetized or baryon-rich jet (for the high 
GeV fluence ratio) from the white dwarf progenitor (for 
the short GeV onset delay) ( ^5^ . 
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